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ABSTRACT 



The presence of non-linear behavior in poten tiodynamic 
polarization plots has resulted in difficulty in determining 
the Tafel constants from such plots, A FORTRAN based program 
involving numerical differentiation techniques using a 
graphical display was used to determine the existence of the 
Tafel regions- 

Varioijs alloys polarized in synthetic seawater and a 
3.5% NaCl solution were analyzed. Although severe 
concentration polarization often dominated the cathodic 
branches the techniques employed did allow for the selection 
of regions which approached linear behavior. The effects of 
concentration polarization in hindering the determination of 
Tafel constants were exemplified by the uncovering of a 
cathodic branch containing a small region where only 
actvivation polarization dominated followed by the onset and 
total domination of concentration polarization. 

A method of determining where the anodic and cathodic 
currents begin to dominate the poten t iodynamic polarization 
curve is introduced. 



3 



TABLE OK CONTENTS 



^r/3-^ 



I, INTRODUCTION 7 

1 1 , DATA TRANSFER 9 

III, CORROSION RATES 13 

IV, GENERAL OVERVIEW OF THE PROGRAM 29 

V, RESULTS AND DISCUSSION 36 

VI. CONCLUSION AND RECOMMENDATIONS 4/ 

LIST OF REFERENCES 49 

APPENDIX A. FIGURES 51 

APPENDIX B. DATA TRANSFER PROGRAM 90 

APPENDIX C, CORROSION PROGRAM 93 

INITIAL DISTRIBUTION LIST 138 



4 



LIST OF SYMBOLS USED 



cm^ 

C.R. 

Ec or- r- 
0CZ r- 

E 

0 

E,W. 

1 cr Of— f— 

id 

i o 

ii_ 

mpy 

mV 

uA 

R 

P- 

P<= 

'7 

7- 

PAR 

PDP 



centimeter squared 
corrosion rate 

corrosion potential of a single metal 
corrosion potential of a single metal 
ox idation potential 
reduc tion potential 
equivalent weight 

corrosion current density of a single metal 
current density 

equilibrium exchange current density 

limiting current density 

mils per year 

mi 1 1 i vol t 

microampere 

Tafel slope 

anodic Tafel slope 

cathodic Tafel slope 

polarization 

acti vation polarization 

concentration polarization 

Princeton Applied Research 

Poten tiodynamic polarization 



5 



ACKNOWLEDGMENT 



Special thanks to my wife Luanne whose patience was 
greater than my perseverance . To Professor Jeff Perkins and 
Mr. Tom Kellogg whose guidance and expertise was most 
appreciated . 



6 



I . iNTEQDyPIlQN 

Recent studies at the Naval Postgraduate School have 
centered on the vibrational charac teristics of certain high- 
damping alloys, [Refs- 1-8]. Most of these studies have 
been concerned with either the measurement of the damping 
capacities of these alloys and/or with the heat treatments 
necessary to attain the microstruc tures associated with the 
observed high— damping capacities. In addition, two separate 
studies were conducted in an attempt to determine the 
corrosion charac ter istics of these high damping alloys in 
marine environmen ts . 

In the studies involving corrosion charac ter istics , both 
laboratory and actual marine environment exposure 
experiments were conducted. The actual marine environment 
exposure tests calculated corrosion rates based on direct 
weight loss methods. The laboratory tests conducted were 
Linear Polarization and Poten tiodynamic Polarization in 
synthetic seawater [Ref. 9] and a 3.5X NaCl solution [Ref. 
10] environments. Although the correlation between the 
laboratory and direct exposure methods was good, the 
individuals involved in both of these laboratory studies 
expressed concern with the difficulty in determining the 
Tafel constants necessary to determine corrosion rates. 
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The actual experiments were conducted on the EG&G 
Princeton Applied Research Model 351 corrosion measurement 
system. The Model 351 utilizes a PAR Model lOOO 
microcomputer featuring two Motorola 60000 microprocessors 
and a touch-screen input to create an environment which is 
well suited for conducting various corrosion tests. The 
system is, however, limited by the lack of a key board and 
more importantly by an operating system which is proprietary 
in nature. These two combine to make any analysis of 
experimental data, other than by the means provided with the 
operating software, impossible unless the data can be 
transfered to another computer. 
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1 1 . DATA TRANSFER 



The PAR 351 and its associated software has the ability 
both to plot the results of an experiment and to send just 
the data points to a printer via an RS 232 serial port. 
Although the software supports only Houston Instruments 
plotters, the printer option is configurable for numerous 
printers as long as the printer supports DCl (XON, ASCII 17) 
and DC3 (XOFF, ASCII 19) protocol. Other RS 232 
requirements such as baud rate, parity, word length, and 
stop bits are configurable as the user sees fit- 

During the initial data transfers, in the present work, 
Hayes Microcomputer Inc. SMARTCOM II modem software was 
used to capture the data. This method was successful most 
of the time but since the two systems were not using modems 
but instead were directly connected via their serial ports 
some problems were experienced and editing of the data file 
after transfer was always necessary- The most serious 
detraction was the need to wade through various menus in 
order to accomplish the transfer. 

In view of the both the problems experienced with using 
SMARTCOM II and the neccessity for the user to be familar 
with the way SMARTCOM II works, it was decided to write a 
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data transfer/capture program which would both speed 
execution and eliminate the need for user knowledge of 
communications so f tware . 

Using the sample TTY data transfer program included in 
the manual for Microsoft GW—BASIC version 2 [Ref. 11] as a 
core a BASIC program was written to accomplish the transfer. 
The core was modified to eliminate the transfer portion and 
most importantly line feeds and carriage returns sent by the 
351 were automatically deleted and a carriage return was 
inserted only at the end of a line containing data. This 
eliminated the problem of blank lines appearing within the 
data file corresponding to page headers and footers. Other 
modifications which were made involved checking for and 
correcting two formatting variations which resulted in 
compatibility problems between the transfered file and its 
use in FORTRAN application programs- The first involved a 
change from a floating format to a exponential format when 
the potential was equal to 0.0000 volts. The 351 normally 
sends potential with 4 significant digits. When the 
potential is 0-0000 volts it is displayed as -XXE-12. The 
formatting shift is not compatible with any FORTRAN formats 
as the total field width changes in addition to the minimum 
number of available digits. This occurance is now trapped 
and the file modified to reflect the actual potential of 
0-0000 volts- The second formatting variation always 
occured and did so in a region which was of extreme interest 
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as follows. When the current density changes from negative 



to positive there is a range where the exponent is E-12. 
Unlike FORTRAN write statements which always reserve a space 
for the - 9 - or - sign, ie. the sign is not printed unless 
specified, the 351 used the space for the - sign if present 
or for the first digit if not. Again this resulted in a 
change in total field width at the occurance of -a.bcdE-12. 
The FORTRAN read statment would read this to be -a.bcdE-1. 
This was handled by simply deleting the least significant 
digit such that the number is read as — a.bcE— 12 when this 
occurs . 

In its final form the data transfer program requires the 
user to reply only a to a (R)eceive or (E)xit prompt and to 
enter a file name if the (R)eceive option is selected. 
Because the 351 can not write a selected configuration to 
disk the user is still required to properly configure the 
351 to the proper settings when powering up. These are: 

PROTOCOL - PRINT 

PARITY - NONE 

BAUD RATE - 2400 

STOP BIT - 1 

WORD LENGTH - 0 

Once these options are properly set on the 351 the 
desired data file should be loaded into the 351 and 
formatted into the linear option. At this point the user 
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selects the (R)eceive option on the Z248, presses the “PRINT 
DATA** prompt on the screen of the 351 and enters a filename 
on the Z248. 

The final version of the data transfer program was 
compiled using Microsoft's BASIC compiler Version l.O. By 
eliminating the writting of the data to the screen in 
addition to the file and using the compiled program vice the 
BASIC interpreter a file which will use 60K of disk space 
can now be transferred in under 5 minutes as opposed to the 
20 minutes it took with SMARTCOM II. A listing of the 
program is included in Appendix B. 



12 



III. CORROSION RATES 



A. THEORETICAL EQUATIONS 

In general the corrosion behavior of a metal can be 
described by the use of the Nernst equation, activation 
potential , concentration potential and mixed potential 
theory . 

The Nernst equation is a modification of the Gibb's 
Free Energy equation and expresses the tendency of a 
corrosion reaction to proceed in terms of EMF based on the 
activities of the products and reactants of the 1/2 cell. 



E = E^-» -* RTloq,(Q) 
nF 



E = Electrode oxidiation potential 
(0 for reduction) 

E*^ = Equilibrium oxidiation potential 
(0 for reduction) 

R = Gas constant 
T = Temperature 

n = Number of electrons taking part 
F = The Farday 

Q = Activities of products/activities of 
reactants 

The potential of the cell is the algebraic sum of 
oxidation and reduction 1/2 cells. 
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Activation polarization describes the current density 



which exists when the cell is polarized, 
imoocj j.€= = ioexp( +ZFaE/RT ) 

ic Atit-«c=>cjxc: “ i c>ex p ( “*ZF ( J. ci)E/RT) 



io = equilibrium exchange current density 
Z = Valance of the active species 

E = Polariztion about the equilibrium potential 
a = fraction of E which applies to i^r^ocjxc 
R,T,F as above 

By manipulating the above equation the activation 
polarization can be written as: 

= 0 logxo (icj/io) 

B = 2.3RT/aZF 

Although activation polarization occurs for both the 
anodic and cathodic reactions the cathodic reaction can also 
be effected by concen tration polarization . When this occurs 
the reduction reaction occurs at such a rate as to decrease 
the concentration of ions in the solution near the surface 
of the metal. This results in a change in the potential of 
the 1/2 cell reaction since the effective activity in the 
Nernst equation changes. 
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The concentration polariztion <= is expressed as: 

^ ” 2 - 3R r 1 og AO (1 “id^ii-) 

nF 

id = operating current density 
±t_ = limiting current density 

The limiting current density is inversely proportional 
to the thickness of the diffusion layer. This thickness is 
dependent on various geometric and environmenta 1 factors and 
can often only be determined experimentally- The effect of 
concentration polarization can be minimized by agitating the 
solution, thus decreasing the layer thickness, or, as in the 
case of oxygen reduction, by heating the solution to reduce 
the concentration of the dissolved gas. Both of these means 
of dealing with concen tration polarization may significantly 
alter the desired corrosion environment and distort the 
correlation between the experimental results and the actual 
in-service corrosion rates. 

In the absence of concen tration polarization a dynamic 
polarization diagram of a cell reaction would appear as in 
Figure 1. The slope of the anodic and cathodic curves 
d0/d(logio id) would be given by B as described in the 
activation polarization expression - 

Most cells are affected by concen tration polarization 
such that the cathodic polarization The 

effect of this is shown in Figure 2- It is obvious from 
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this that log^o (id) approaches logxo (it,) the linear 
behavior of the polarization curve is distorted and in cases 
where it, and io are not sufficiently separate, a linear 
region may not exist- 

Mixed potential theory concerns concerns the immersion 
of a single piece of metal into an electrolyte and the 
subsequent reactions. The theory consists of two simple 
hypotheses . 

1- Any electrochemical reaction can be divit’^ed into two or 
more partial oxidiation and reduction reactions. 

2. There can be no net accumulation of electrical charge 
during an electrochemical reaction. [Ref. 12 : pp 314] 

It is these hypotheses that allows for the 
determination of two important parameters necessary for 
estimating corrosion rates. When the cathodic branch of the 
1/2 cell with the highest equilibrium potential intersects 
the anodic branch of the other 1/2 cell the corrosion 
potential 0oc*r-r- and the corrosion current density icor-i- can 
be used to estimate corrosion rates. 

When the corrosion system consists of only a single 
metal in contact with a solution which contains a single 
reduction/oxidation reaction the overall reaction is the 
linear sum of the reduction currents and the oxidation 
curren ts . 
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If the solution contains more than one 



reduction/oxidation reaction the overall reaction may 
proceed as the linear combination provided that the 
equilibrium exchange current densities of the secondary 
reduction reactions are within one order of magnitude of the 
main reduction reaction. The secondary oxidiation reactions 
proceed in the same manner. As a result the secondary redox 
reactions affect either, both, or none of the branches of 
the overall reduction/oxidiation reaction. [Ref. 13: pp. 

199] 

If concentration polarization is not present both the 
anodic and cathodic polarization curves would consist of a 
of horizon tial lines, representing the current added by the 
reaction at that potential, and slanted lines of some slope 
8 representing the activation polarization of all the 
oxidation or reduction reactions which are influential at 
that potential. In the presence of concentrat ion 
polarization the possibility of introducing numerous non- 
linearities into the cathodic branch is distinct. As the 
overall reduction reaction is the sum of all the influential 
reduction currents any region could consist of zones where 
the activation polarization of the ith reaction dominates 
the concentration polarization of the (i— l)th reaction. 

This can make the poten tiodynamic curve extremely difficult 
to analyze. 
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B. DETERMINATION OF CORROSION RATES 



As previously described the presence of numerous 
electrodes and/or metals in the corrosion system can 
significantly alter the shape of the polarization curve. 
Despite this if it was possible to design a current 
measurement system that was capable of measuring only the 
anodic and cathodic currents the determina t ion of i^or-r- 
could still be simply made since log^o i<=oi-r- would be at the 
intersection of the anodic and cathodic branches of the 
polarization diagram. 

In the real world polarization curves take on a 'Y' 
shape as a result measuring instrument limitations and the 
complexity of the corrosion environment. 0c=oi-i- can be 
measured directly since it is the equilibrium potential of 
the system. As the potential is scanned towards 0c=c3i-r- the 
net current (i«r»ocJAc “ ic«^r^oa xc ) is measured. Only when 
the potential is displaced far enough from does one of 

the two currents overwhelm the other and does typical Tafel 
behavior take place. 

The presence of numerous 1/2 cell reactions, 
concentration polarization and the need for one of the 
currents to overwhelm the other can all combine to make the 
determina tion of ice^»-f- difficult. 

If concentration polarization is not present at some 60 
from 0CO.-I- the polarization diagram should adopt true linear 
behavior. The amount of 80 at which the curve should adopt 
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linear behavior is in itself open to disscussion . Ailor 
states that S0 should be at least 50mV from 0c<=,r-r- CRef. 14: 
pp. 199]. Once the linear regions are found two methods can 
be employed to determine i<=c=»f-.- - The simplest is to draw two 
straight lines from the linear regions towards 0crc3i-r-. The 
two lines should intersect at ic=of-t-- A second method is to 
measure the slopes B (Tafel constants) of the linear regions 
of the anodic and cathodic branches. These values can be 
used in conjunction with data determined from a polarization 
resistance test to determine the value of icor-r- by use of 
the equation: 

i.= c=,r-r- = 

2.3(B^+Bc=)PR ' 

f 

If concentration polarization is predominant over the 
cathodic region of the polarization curve, the so-called 
'knee' method is often employed- In this method a tangent is 
drawn on the knee of the cathodic branch such that the 
angles formed on both sides of the tangent with the cathodic 
branch are equal. When this has been accomplished either of 
the methods discussed in the previous paragraph may be used. 

Many other methods have been developed to determine 
ic=c 3 r-r- or to determine corrosion rates directly from 
polarization data. Princeton Applied Research uses a Chi— 
squared minimization technique. This technique is included 
in the software package which operates their Model 351 
corrosion measurement system as PARCalc. By by their own 
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admission, since the program is terminated based on a change 
in the average value of the sguare of the relative 
deviations between two successive iterations, it is 
necessary to determine the Tafel region before starting the 
process. [Ref. 15: pp IX-10] 

LeRoy [Ref. 16: pp 1006-10121 postulated the 
determination of Tafel Slopes based on polarization 
resistance technigues. Without repeating the entire 
discussion it is sufficient to say his method of determining 
Tafel slopes based on expected values from activation and 
concen tration polarization expressions, based on the 
assumption that the parameters necessary to evaluate those 
expressions were available, did not accurately describe the 
observed values. [Ref- 17: pp 1988—1989] 

It would also seem possible to solve the overvoltage 
expression for the terms in guestion. Recalling the 
overvoltage expression : 

0 - 0c=or-r- = Be 1 og xo ( ici / i c ) 2.3RT loqxo( 1 - i<d/it.) 

nF 

The problem is that Bt=, io and it_ all represent unknowns 
and the ix. term is located in such a manner as to complicate 
the solution. 
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The logio (i - i<a/iL_) term can be replaced by the series 
expansion : 

log^od + X) = 2.3C X - + X=" - X^ ... 

2 3 4 n 

where X = “i<a/ii» 

The final equation to be solved would be: 

0 - 0^«3r-^ = Bc=logxo(ic)-Bc:logxo(ic)+2.3C2(~1)--^.X-] 

n 

The difficulty with this approach is that a 
sufficient number of terms must be carried in the series in 
order to accurately represent the logxo term. This causes 
two severe problems. The first is n+2 data points must be 
used to solve the system of equations. When one considers 
that id is normally on the order of ^/cm^, exponentiating 
it to the nth degree can easily result in a matrix to be 
solved involving coefficients of 50 or 60 orders of 
magnitude. The end result is that while the system can be 
solved such that the coefficients determined will result in 
the accurate represen tation of the actual overvoltage, the 
coefficients themselves have been influenced by round-off 
error and catastropic cance 1 lation , and represent merely 
numbers. 

The other problem exists with attempting a solution of 
this form. As discussed previously with the possibility of 
mixed electrodes, the associated redox processes may or may 
not effect the overall behavior of the system. Even if the 
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limiting current densities ii_ were all known, it would still 
be necessary to determine their effect on the overall 
process in order to include them in the system of equations 
to be solved. As a result, for m different electrodes, m 
experiments involving just the electrode and the metal would 
have to be carried out just to determine their influence. 
Consider the use of such a method with synthetic seawater. 

At least 10 different chemical substances could effect the 
overall reaction in addition to the major metal oxidiation 
and hydrogen reduction reactions. It should be obvious that 
if the limiting current densities are not known the size of 
the system of equations to be solved could easily exceed the 
number of available data points. More importantly, one 
would have to assess the advantages of performing a 
polarization in a mixed medium if the same test has to be 
performed separately on the individual components first- 

C. CALCULATION OF TAFEL CONSTANTS 

In order to determine the location, if one existed, of 
the anodic and cathodic Tafel regions two numerical 
dif f eren tation techniques were employed. These two methods 
were the Four Point Central Difference method and the use of 
a Cubic Spline interpolating polynomial. A graphical 
display of the derivatives allowed the user to determine the 
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region of linearity. Once this region was determined a 
first order linear regression is performed on the original 
data within the region to determine the Tafel slopes, 
i - Four P oint Ce ntra l Di fference Method 

When the data points are egually spaced the Central 
difference method may be used to approximate the derivatives 
at a point. The use of this method can result in 
considerable time savings over the cubic spline method as 
the derivative can be calculated without having to resort to 
the solving of a system of equations for each point in 
question. The method is not as accurate as the cubic spline 
method and graphical represen ta tions of the derivatives will 
often contain more "noise” than those derived from the cubic 
spline method. 

The four point Central Difference formula for 
computing derivatives is: 



d0 


= 12h 




d(logio id) 


8( fx - f_x) + 


( f-:z - ) 


Where: 


h = 


0± ~ 




fr, = 


logxo(ioi) at the 
1 Og iO ( icai ) I 


Ith location from 



The Central Difference method has an error term 
associated with it. In this case the error term is on the 
order of 1/h^. Although this may appear to be a substantial 
term, it should be remembered it is a possible error, the 
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method is only being used to determine the Tafel region and 
when used the Tafel constant is inverted so the error term 
would actually be h^ . 

2 - Cubic S pl ine Inte rpolating Po l yno mial 

When trying to numerically determine derivatives of 
experimental data which may or may not be equally spaced the 
use of interpolating polynomials is often suggusted. One of 
the most accurate methods is to pass a cubic spline thru the 
data points and then to dif f eren tiate the resulting equation 
to arrive at the derivative at a desired point. CRef. 18: 

PP- 242] 

The advantage of using a cubic spline is that in 
deriving the system of equations to arrive at the constants 
a, b, c, and d which form the equation: 

y = aj(x - Xi)^-^- bi(x ~ Xi)^ + Ct(x - xi) + dx 
the system is formed in such a manner that the value 
of the function, first derivative and second derivative are 
the same for the pair of cubics which join at each point - 
This requires that: 

y = 3ai ( X — Xj)^ 2bx(x — Xi) + Cx 
and y' ' = 6ai ( x ~ Xj) + 2bi 

In order to simplify the the mathematics involved 
the equations for the system are written in terms of the 
second derivatives of the in terpolating cubic St . 
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The resulting interpolating cubic will then take on 
the form of : 

h I — J.S I — i+2S X ( h I — X ) "^h I S X ( y X *♦* i ) ~ ( y i ~y I — j. ) 

hi hx— i 

where hx = x - Xi 



Solving the system leads to the following operations 
to determine the values of ax , bi , Ci and c/x - 

aj — (Si-^x ■“ Sx)/^hi 
bx = Sx/2 

Cz = ((yi^j. - yi)/hi) - ((2hiSi + hxSx^i)/6) 

dj. = yi 



The resulting system of n - 2 equations in Si 
involves n pairs of data points in order to generate the 
required number of equations. To arrive at the two 
additional equations to solve for Sj. and Sr», constraints are 
specified which pertain to the conditions at the ends of the 
curves. The three choices for the end conditions are: 



1. Si = Sr, = O, which implies that the end cubics 
approach linearity at their extremities. 

2. Si = Ss?, Sr, = Sr,-i- This assumes that the end cubics 
approach parabolas at their extremities. 

3- Si is a linear extrapolation of Si and Sz- Sr, is a 
linear extrapolation of Sr , -2 and Sr,-i - 
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In this case since the interest is in determining 
the derivatives d0/d(logjto icj ) pointwise, and not in 
developing an interpolating polynomial, it was decided to 
use the end conditions of Sx = = O, and to use enough 

data points as to force the inconsistences caused by the 
assumed end conditions away from the point in question. 

Four data points on each side of the point in question were 
used in determining the the derivative at the point in 
question. This resulted in the solving of a 7 x 7 system of 
equations for each data point. By locating the point in 
question in the center of the data set, the accuracy of the 
estimate involved was improved, and the mathematics involved 
in determining the derivative were simply solving for c^. . 

^ • Linear Re g ress ion 

Once the anodic and cathodic Tafel regions have been 
determined, the data points which comprise these regions are 
used in a least squares curve fitting routine to determine 
an equation of the form: 

0 = ^(logxo (id)) 

where P represents the Tafel slope. The values of a 
and b are found by solving the two following equations 
simu 1 taneous 1 y . 

( 1 og jto ( i d ) I ) ^ ( 1 og xo ( i d ) I ) “ 2( logxo( id ) i ) (0i ) 

( 1 og xo ( id ) X ) “*■ = Z( 0 i) 

where N is the number of points in the Tafel region. 
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When the PAR 351 Corrosion Measurement system scans 
the potential during a Poten tiodynamic test the current is 
measured at intervals of 2mV or 0*5mV/ of potential, 
depending on the relative change in the current. The 
sampling rate is normally 0.5mV when the current shows a 
substantial change and 2mV when it does not- This usually 
results in a sampling rate of O.SmV in the Tafel regions, 
although severe concen tration polarization in the vicinity 
of icrof-r- may result in a limited number of points having 
0.5mV intervals on the cathodic branch. Two other aspects 
may also effect the potential intervals at which the current 
appears to have been sampled. 

When data is sent out the serial port, the 351 
normally uses a format similar to the FORTRAN Ga.b type. 

Not counting leading zeros, four significant figures are 
used. This results in a problem in the recorded potential 
values- When the potential is greater than or equal to i.O 
volts, or less than or equal to —1-0 volts, only three 
decimal places are carried. As a result, when the values of 
fall in this range and the scan rate is O.SmO, the 
recorded data seems to reflect that the potential has not 
changed for two successive current measuremen ts - 

Since the potential is being scanned it is rather 
obvious when this occurs that the recorded data does not 
accurately represent the experimental data. Since this 
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pattern is easily recognized, the array containing the data 
can be manipulated to more accurately reflect the 
exper imen ta I data • 

Occasionally the potentials recorded display no 
regular pattern and the data, although questionable in 
nature, must be used as recorded. 

When calculating the derivatives, the anodic and 
cathodic branches are performed separately- The 
calculations begin at the third data point away from Ecoi-r- 
when the Central Difference method is used, and at the fifth 
when using the cubic spline method. The calculations 
continue as long as the data points continue to suggest a 
slope of of the expected sign (+/—). This automatical ly 
excludes the passivating region of the anodic section, as xt 
should since it is not part of the Tafel region, and 
excludes the possibility of the program terminating 
execution if it were to calculate an infinite slope. In 
general this strategy worked as expected and provided enough 
data points to properly examine the Tafel regions. The only 
experiment in which it failed was for the SONOSTON sample 
tested in 3-5Z NaCl solution. As can be seen by examining 
the Potentiodynamic curve of this sample (Figure 3), the 
cathodic branch exhibited behavior throughout which caused 
either method to terminate execution without generating any 
derivatives- The reason for this behavior is unknown. 
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IV- GENERAL OVERVIEW OF THE PROGRAM 



The program was developed so as to allow for its use 
without the supposition that the user was familar with the 
FORTRAN langauge or with personal computers in general. 
Numerous efforts were made to minimze the possibilty of the 
program prematurely terminating execution in the event of 
improper input or through floating-point operations 
resulting in invalid expressions. By taking advantage of 
the color capability of the monitor screen, colors are 
changed with respect to the type of input expected, and to 
display messages resulting from either input errors, data 
file format errors or a failure to properly follow the steps 
necessary for proper program execution. 

The program was written using Microsoft FORTRAN 
Optimizing Compiler version 4.0. Graphics were achieved 
using Microcompatibles, Inc.'s GRAFMATIC library. Plotter 
support was achieved using Microcompatibles, Inc.'s 
PLOTMATIC library. 

A summary of the main program and its subroutines is 
given below. The main program and its called subroutines 
are listed in Appendix C. The subroutines are presented in 
the order in which they appear in the main program- 
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A. START -FOR 

This is the main program. The name of the first data 
file to be used is entered. Following calls to DATAIN, 

CHECK and DATAM a menu consisting of 6 options is presented: 

1. Plot the poten tiodynamic curve 

2- Generate Tafel Slopes using cubic spline 

3. Generate Tafel slopes using central difference 

4. Use another data file 

5- Overlay two potiodynamic curves 

6. Exit 

Option 1 generates calls to GRAPHl and PLOT 

Option 2 calls SLOPEl . Following this a first orde»^ 
linear regression of the anodic and cathodic branches using 
locations within the original array determined during a call 
to GRAPHl is performed and the results along with original 
poten tiodynamic curve are passed to GRAPHl - 

Option 3 calls SLOPE- Upon return from SLOPE this 
option proceeds as Option 2 

Option 4 returns the user to the 3rd executable line of 
START. 

Option 5 calls DATAIN, CHECK and DATAM while retaining 
the data from the first file. 

Option 6 obviously exits the program. 

B. DATA IN. FOR 

This opens the input file and, using a formatted read 
statement, generates the array containing potential and 
current density. Additionally this subroutine counts the 
total number of lines in the file. This is done to ensure 
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that formatting changes within the data file which were not 



corrected during the data file transfer will not have not 
caused the array generation to prematurely terminate- If a 
formatting error has been detected, an error message is 
displayed in red and the line number within the file is 
identified. This allows the user to exit the program and 
correct the error before continuing. 

C. DAT AM. FOR 

This converts the absolute value of the current density 
to its loqxo value and determines the position of Ec=or-r- with 
the array. 

D. CHECK. FOR 

This subroutine determines the number of data points in 
the anodic and cathodic branches to be used in the 
determining the derivatives when using the central 
difference method. This is done by checking log current 
density for values which might result in the calculation of 
an infinite slope (normally the onset of pitting or 
passivation) and for equal increments in the recorded 
potentials. Occasionally the recorded potentials will not 
reflect the polarization which has taken place. This occurs 
when the absolute value of the recorded potential is equal 
to or greater than 1.0. In these ranges the O.SmV voltage 
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increment falls outside of the total width of the field and 



is not recorded. As this is a recognizable pattern, at the 
users option the problem can be corrected. 

E. SLOPE. FOR 

This is the first subroutine called when using the 
Central Difference method. It assigns graph titles, axes 
labels and calls DATADEL-FOR. 

F. DATADEL-FOR 

Using the points determined in CHECK, this calculates 
the numerical derivatives using the Central Difference 
method . 

G. SLOPE 1. FOR 

This assigns graph titles and labels. It also 
determines the number of points, MI , to be used on both 
branches in the cubic spline method by checking only for 
values of logxo current density which could result in the 
calculation of slopes of infinite value. It calls CSPLIN. 
separately for the anodic and cathodic branches. MI anodic 
is not necessarily egual to MI cathodic. 

H. CSPLIN. FOR 

For MI number of points, this generates the 7x7 matrix 
and the 7x1 vector of values described in the cubic spline 
section MI-‘9 times. As each matrix is formed it is solved 
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by a call to LINSYi.FOR. Upon return from LINSYl the 7x1 
vector contains the solution which is used to calculate the 
deri vative. 

I. LINSYI.FOR 

Using partial pivoting this subroutine solves the 
augmented 7x8 array by Gaussian Elimination. The original 
7x1 vector contains the solution [Ref. 19]. 

J. GRAPHl.FOR 

This is the graphics output of the selected main menu 
option. Through the use of the numeric keypad as softkeys 
the user can scale the axes as desired. The option to 
return to the original graph or the previous display always 
exists. When plotting the results of either numeric 
dif ferentation method, only those points which consecutively 
fall within the desired region are plotted. Once the axes 
have been scaled such that graph contains only the linear 
region of the anodic or cathodic branches the position of 
the start and stop points within the original array 
containing the derivatives are saved. Since, when 
calculating the derivatives, the data point containing Ecorr 
is always used as a reference, ie. (depending on the method) 
the first derivative for each initial branch is either 3 or 
5 points away from Ecc»r-r-» the start and stop points also 
represent the position of the Tafel regions within the 
original potential and loglO current density array. When 
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the graph is not of the numerical derivatives, any data 
point which falls within the desired axes is plotted. If 
desired GRAPHl calls PLOT. FOR for a hard copy output of what 
is currently displayed on the screen. 

K. PLOT. FOR 

This plots the graph displayed on the screen. Although 
many of the arguments passed from 6RAPH1 are used directly 
as passed, this is not a 'screen dump'. The resolution of 
the plctLjr is much higher than that of the screen. As a 
result what appears to be a one pixel Jump on the screen 
will not plot as such. Two sets of tabulated data may be 
plotted on each set of axes. At the users discretion 
different color pens may be used to plot the axes, first 
output, second output, and labels. 

L. CONFIGURATION 

In its final form, the program has been compilied and 
linked into an executable file which runs on a Zenith 240 
computer which is equipped with a 80286 microprocessor , an 
80287 Numeric Processor Extension, a EGA graphics card with 
64K video RAM, and a color monitor which supports the EGA 
color mode. The plotter used is an Houston Instruments DMP- 
40 series digital plotter. 

The plotter library is plotter series specific. As a 
result of this, the use of a plotter other than a Houston 
Instrument's DMP— 29 or greater would require obtaining 
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another library from Microcompatibles Inc. As many of the 
subroutines contained within PLOTMATICS are similar to those 
called within GRAFMATICS, it can probably be safely assumed 
that little if any editing of the subroutine PLOT. FOR would 
be required. 

To use the program with a video configuration other than 
that previously described would require editing of the 
subroutines PLOT. FOR, GRAPHl.FOR, DATAIN.FOR, START-FOR and 
CHECK. FOR. DATA IN, START and CHECK would require editing 
only if the video configuration was monochrome- PLOT and 
GRAPHl would require editing for any graphics configuration 
other than EGA. 
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V 



RESULTS AND DISCUSSION 



A. SYNTHETIC SEAWATER 

Table 1 below compares the Tafel constants determined by 
using the new program with those obtained by Escue [Ref. 9]. 



TABLE 1 

TAFEL CONSTANTS OBTAINED FROM 
SYNTHETIC SEAWATER EXPOSURE 



MATERIAL 


0. 


0-* 


0c= 


0=* 


Ti-50Ni 


0.1400 


0.174 


0.1155 


0.132 


1020 Steel 


0.0758 


0.068 


0.1481 


0.114 


304 Steel 


0.0241 


0.178 


0.0207 




7075 A1 


0.3790 


0.087 


0.1584 


0.062 


Fe-Cr-Mo 


0.1775 


0.121 


0.1226 


0.141 


Fe-Cr-Al 


0.1840 


0.180 


0.1369 


0.151 


Cu— Mn— Al— Fe— 


Ni 0.0399 


0.048 


0.1579 


0.120 


Cu-Mn-Al 


0.0633 


0.035 


0.0652 


0.045 


Cu-Zn-Al 


O . 0804 


0.062 


0.0785 


0.075 


630 Bronze 
» From Escue 


0.0403 
[Ref. 9] 


0.031 


0.0724 


0.041 



Table 2 below compares the current densities and 
related corrosion rates of the samples tested in synthetic 
seawater- The first line for each material represents the 
current density as determined by the intersection of the 
Tafel lines- In many of the cases the anodic and cathodic 
Tafel lines did not intersect at a single point and are 
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joined by a horizontal line. In these cases icor-r* was 
determined by taking the average of the anodic and cathodic 
icor-r- at the point where the lines intersected a horxzontal 
line form by 0c=or-r-- 



The second line represents the results obtained by 



combining the Tafel slopes with the polarization resistance 



da ta • 



TABLE 2 



CURRENT DENSITY AND CORROSION RATES 



MATERIAL 




i <r cr> r- r- * 


CR 


CR* 




uA/cm^ 


uA/cm^ 


mpy 


mpy 


Ti - 507.N1 


1.245 


1.734 


0.573 


0.799 




4.182 


4.960 


1.822 


2.160 


1020 Steel 


4.212 


3.837 


1.981 


1.805 




20.325 


17.201 


9.580 


8.110 


7075 A1 


1 .413 


0.421 


0.591 


0.176 




3.067 


0.994 


1.281 


0.415 


304 Steel 


0.106 


0.728 


0.048 


0.329 




0.157 


1 . 190 


0.071 


0.537 


Fe-Cr — Mo 


0.750 


0.804 


0.354 


0.380 




2.294 


2.060 


1.085 


0.974 


Fe-Cr-Al 


0.822 


0.830 


0.362 


0.365 




1.183 


1.240 


0.520 


0. 545 


Cu— Mn— A 1 — Fe— Ni 


8.241 


8.854 


4.165 


4.070 




9.662 


10.400 


4.882 


5.260 


Cu-Mn-Al 


2.089 


1.119 


1.111 


0.595 




1.459 


0.894 


0.784 


0.481c 


Cu— Zn-Al 


2.239 


1.897 


1.136 


0.963 




3.780 


3.230 


1.919 


1.640 


630 Bronze 


1.762 


1.442 


0.901 


0.737 




6.878 


4.690 


3.550 


2.420 



* From Escue [Ref. 9] 
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i . T i - 50y.N (TIN I LG ) 

These results compared favorably with those obtained 
by Escue. Figure 4 shows the derivatives as calculated by 
the cubic spline method. As can be seen in this figure no 
true linear regions are distinguishable in either of the two 
branches- When this occurs the "knee” method can be 
employed by examining the branches for a range where the 
derivatives oscillate about a linear midpoint. In this case 
both the anodic and cathodic branches showed this 
oscillation in the log current density range of -5.20 to - 
5.51. Figure 5 shows the resultant Tafel lines plotted on 
the original poten tiodynamic curve. Although the curves 
intersect 0c:cDr-r- at what could be considered to be a single 
point the icor-.- determined in this manner differs greatly 
from that calculated by using the Tafel constants and the 
polarization resistance data. 

2. iPJZ^Carb^ 

Figures 6 and 7 show the results as obtained for the 
1020 carbon steel sample. Of all of the samples tested this 
one exhibited a cathodic branch with the least tendency to 
behave in a linear fashion. If there was any linear 
behavior in the cathodic branch it occured at a log^o iccD*-*- 
of about -5.2. The anodic branch displayed a much longer 
lasting region around logi.o icror-r- = ~ 4.3. This sample also 
showed the greatest disparity in the value of ic<=»r-f- 
depending on the manner in which it was calculated. 



38 



3 . 






In this case the central difference method showed a 
linear region in both the anodic and cathodic branches. As 
can be seen in Figure 8 the linear regions exist at the 
pronounced peaks in the derivatives- In this case the Tatel 
lines did intersect at a single point as seen in Figure 9. 

^ • 304 Stainless Steel ( SST30 4LG) 

The calculated corrosion rates and the respective 
Tafel constants differed greatly from those obtained by 
Escue- Although the calculated corrosion rates correlated 
better to the actual weight loss methods, the need to 
terminate the calculating of derivatives before the onset of 
passivity may have caused a linear region in the anodic 
branch to be omitted. These points can be seen by the sharp 
peak on the anodic branch on Figure 10, and the 
corresponding point on the poten tiodynamic curve Figure 11. 
The Tafel lines intersected at a single value of logi.o icor-r- 
but this time there was a good correlation between the 
calculated and graphical icor-r-* 

5 - F e-Cr- Mo (VACRQILG) 

Again the presence of any truly linear regions is 
questionable, although both branches. Figure 12, show a 
range where the deviations from linear behavior are minimal. 
Figure 13 once again shows a theoretical intersection of the 
Tafel lines. 
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^ LyACR02L^ ) 

This alloy showed behavior similar to the Fe-Cr-Mo 
alloy described above. What may be considered a truly 
linear region in the cathodic branch is shown in Figure 14 
just to the left of the sharp peak. The overlay of the 
Tafel lines on the potentiodynamic curve is shown in Figure 
15- The Tafel constants derived show a closer relationship 
to those calculated by Escue, than did a comparision of the 
Fe— Cr~Mo sample. 

^ - Mnr A 1 -Fe-Ni (S0N0STL6) 

SONOSTON, Figures 16 and 17, displayed a linear 
anodic branch and, as did most of the alloys, a cathodic 
branch effected by concen tra tion po 1 ar ization . 

8 . Cu--rin -A 1 _ 

The result here again showed a linear anodic branch 
but in this case the concentration polarization did not 
overwhelm the cathodic branch in such a manner as to distort 
the region where activation polarization is predominant- 
Figure 18 shows the graph of all the derivatives and Figures 
19 and 20 show the graphs of the derivatives in what was 
decided to be the Tafel region. Figure 21 shows the 
original potentiodynamic and the Tafel lines. 

9 . Cut Z n -_A 1__ jJDLAL^^ 

Experinced similar behavior as the Cu— Mn— A1 sample, 
□nee again a linear region was very apparent in the anodic 
branch. Figures 22 and 23. 
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i O - 630 Series Bronze (BRNZL6 ) 



This was the only sample of those tested in the 
synthetic seawater where the Tafel regions for the anodic 
and cathodic branches were found in the same region of log 
icor-r-. As seen in Figure 24 concentration polarization 
again quickly dominates the cathodic branch. Figure 25 
displays the overlay- 
B- 3.57. NaCl SOLUTION 

In a similar fashion Table 3 lists the Tafel constants 
and Table 4 the resultant current densities obtained from 
the samples tested in the 3.57 NaCl solution. In general 
the current results show a stronger correlation with the 
corrosion rates obtained by Akthar [Ref. 10] than did the 
comparision of the current results with those obtained by 
Escue [Ref. 9] for the synthetic seawater. 

TABLE 3 

TAFEL CONSTANTS OBTAINED FROM 



3. 


5 7. NaCl 


SOLUTION 


EXPOSURE 




MATERIAL 


0. 


B-* 


Be 


Bc=» 


304 Steel 


0.1259 


0.4117 


0.0951 


0.117 


7075 A1 


0.0116 


0.010 


0.0383 


2.36 


Fe-Cr — Mo 


0.2821 


0.234 


0.1053 


0.110 


Fe-Cr-Al 


0.1203 


0.214 


0.1292 


0.1285 


Cu-Mn-Al 


0.0082 


0.0137 


0.0337 


1.379 


630 Bronze 
♦ From Akthar 


0.0524 
CRef. 103 


0.570 


0.1472 


0.321 
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TABLE 4 



CURRENT DENSITY AND CORROSION RATES 



MATERIAL 


■l-cor-r- 


icor-r- * 


CR 


CR* 




uA/cm= 


uA/cm= 


mpy 


mpy 


7075 A1 


6.960 


7.406 


0.048 


2.965 




10.460 


11.700 


1.842 


4.685 


304 Steel 


0.052 


0.100 


0.026 


0.050 




0.012 


0.020 


0.006 


0.010 


Fe— Cr — Mo 


0.043 


0.040 


0.019 


0.018 




0.609 


0.613 


0.274 


0.276 


Fe-Cr-Al 


0.493 


0.660 


0.234 


0.313 




0.916 


1.180 


0.428 


0.551 


Cu-Mn-Al 


2.864 


6.062 


1.454 


3.078 




5.251 


10.800 


2.666 


5.483 


630 Bronze 


0.093 


0.129 


0.048 


0.066 




3.570 


4.480 


1.842 


2.312 



♦From Akthar QRef. 10] 

1 - 707 5 Alu minum (AL7075) 

In this case the effect of concen tra tion 
polarization made any estimate of the cathodic Tafel 
constant extremely difficult. As seen in Figures 26 and 27 
concentration polarization completely overwhelmed the 
cathodic branch- As a result only the "knee** method cuold 
be used, but its accuarcy is questionable. The anodic 
branch displayed typical linear behavior. 

2- 304 St ainle ss S teel (S ST304) 

The cathodic branch of this sample did display a 
linear region located between the sharp peaks in Figure 28. 
The anodic branch displayed numerous attempts to passivify 
as seen in Figure 29. 
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3. 



Fe-Cr~M o (VCMOPD) 



Similar behavior to the 304 Stainless steel above. 
The bandwidth exhibited by the cathodic branch near log^o 
icot-f- of about - 6.5 Figure 30, typifies the knee in the 
curve. Again numerous attempts to passivity limited the 
range of data which could be use in the anodic branch Figure 
31. 

4. Fe-Cr-Al ( VACAL2 ) 

Both the anodic and cathodic regions. Figure 32, 
displayed a region where the behavior could be considered to 
be linear. Although the Fe— CR— Mo and Fe— Cr— A1 samples 
tested in synthetic seawater displayed similar behavior, the 
Fe-Cr — A1 sample tested in the 3.5% NaCl solution did not 
exhibit the tendency to passivity. Figure 33, as did the Fe- 
Cr — Mo sample. 

5. Cu-Mn -A 1 (1 NCRMTE ) 

This sample best typifies the problem associated 
with concentration polarization. On Figure 34 a linear 
region can be seen in the cathodic branch at a log^o icoi-r- 
of about -5.3. This region is followed by a rapid decrease 
in the slope of the curve indicative of concentration 
polarization. When the original poten tiodynamic curve is 
viewed thru out its entire polarization range. Figure 35, 
very little can be determined due to the great changes in 
current density on the anodic branch. When the range in the 
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cathodic branch displaying the linear bevaior is scaled 
properly a textbook example of concentration polarization is 
seen. Figure 36- The overlay is shown in Figure 37. 

^ • 63 0 Series Bronze (BRZ630PD ) 

As seen on Figures 38 and 39 this sample experienced 
similar behavior to the one above- All though the linear 
region directly before concentration polarization is not as 
well defined the cathodic behavior is dominated first by 
activation polarization and then by concentration 
po lar i za tion . 

C- COMPARISION OF CURRENT DENSITIES 

In many of the samples tested the corrosion rates as 
determined by the intersection of the Tafel lines were 
significantly lower than those obtained by the use of the 
Tafel Slopes and linear polarization data. As previously 
stated the Tafel regions should begin to dominant the 
polarization curve at an polarization increment of about 
50mV from Ec=or-r-- This implies that a 1 inear region should 
not be present at increments less than this. With these 
tests a linear region was often found within this zone. 
Because of the n^eneral shape of a poten tiodynamic 
polarization curve the resultant lines drawn from this 
region will always have a lower slope and thus intersect at 
a lower value of i. The disparity underscores the need to 
use the average of icor-r- as determined by the intersection 
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of the Tafel lines and that calculated by the Tafel slopes 



and linear polarization data when estimating corrosion 
ra tes . 

Another interesting feature was noticed when conducting 
some data manipulation of the potentiodynamic polariza tions 
analyzed. When the absolute value of (0 - 0c=or-r-) was 
plotted against logj.o icj in many of the cases the anodic and 
cathodic branches began to deviate significantly at a log id 
corresponding to the intersection of the two Tafel lines. 
This may be coinicidental but also may warrant further 
discussion. When a linear polarization resistance test is 
conducted the potential is scanned with a range of 0c=cDr-r- +/“ 
25 mV. The polarization resistance is determined by 
calculating the slope of the linear region about i = O 
A/cm2. If a linear polarization curve is plotted in the 
same manner as above, while neglecting the singularity at i 
= O A/cm^, the deviation discussed would determine the end 
of the linear region. Returning to the equation used to 
determine icor-r- using the Tafel constants and the 
polarization resistance it can be shown that in order to the 
polarizatxon resistance to be a constant that 0*0*= = 
C<(0*+0c=), where C is a constant of proportional i ty . The 
current being measured during a potentiodynamic polarization 
is equal to the difference between the anodic and cathodic 
currents. Only when one of the currents begins to overwhelm 
the other would there be a substantial change in the 
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measured current with respect to a given overvoltage- This 
would correlate to the deviation between the anodic and 
cathodic branches on a plot of the absolute value of (0 — 
0c=or'r*) vs 1 og ( icj ) • It also should correspond to the end 
of the linear region on the linear polarization curve since 
at this point the constant of proportional i ty would no 
longer be a constant as one of the two branches is beginning 
to overwhelm the other- It is postulated that a technique 
of this method could be used to determine the minumium 
overvoltage at which the search for a Tafel constant could 
begin. It is recognized that this does not account for the 
possibility of completely symetric anodic and cathodic 
branches but if this was the case there would be no problem 
in resorting to the 50mV rule proposed by Ailor- 
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V I • CONCLUS JLQN AND RECpMMENpAT I ONS 

The FORTRAN based graphics program proved to be 
extremely valuable in determining the existence of linear 
Tafel regions particularly in those cases where 
concentration polarization dominated the cathodic branch, 

With the ability to transfer the experimental data to a 
programable computer the opportunity exists to further 
analyze the data. This could allow for the determination of 
which electrodes are effecting the main redox reaction when 
the environment consists of a solution of mixed electrodes. 
In order to accomplish this it is recommended that a single 
high damping alloy be subject to individual po ten tiodynamic 
polar izations containing a single component, at its normal 
concen tration , of synthetic seawater- Through the use of 
graphical techniques it should be possible to isolate the 
influential redox reactions- 

In order to accomplish this it will be necessary to 
conduct all the experiments in solutions which contain the 
same amount, or at least a monitored amount, of dissolved 
oxygen. This is one factor which was not considered in the 
previous two experiments and may have made the analysis more 
difficult than necessary. 
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In lieu of this approach an effort should be 
determine the limiting current densities causing 
concentration polarization . 



made to 
the severe 
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Figure 1. Theoretical activation polarization curve 
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Figure 2. Theoretical activation and concentration polarization curve. 
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LOG I R/cm2 

Figure 3. PDP plot of Cu-Mn-Al-Fe-Ni alloy in 3.5X NaCl solution. 
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Figure 4. Cubic Spline derivatives of Ti-50% Ni alloy in synthetic 
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Figure 5. Taiel lines overlaid on PDF plot of Ti-507. 
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POTENTIODYNflMIC 51Q20 lG 
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POTENTIODYNflMIC 55T304LG 
TflFEL SLOPES: CENTRAL DIFFERENCE 
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Figure 10. Central Difference derivatives of 304 Stainless steel in 
synthetic seawater . 
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Figure 11. Tafel lines overlaid on PDP plot of 304 Stainless steel in 
synthetic seawater. 
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POTENTIODYNRMIC VRCROILG 
Bc= 00.1225 Ba= 00.1775 
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POTENTIODYNRMIC VRCR02LG 
TflFEL SLOPES; CUBIC SPLINE 
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Figure 14. Cubic Spline derivatives of Fe-Cr-Al alloy in synthetic 
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POTENT lODYNRMIC 50N05TLG 
TRFEL SLOPES; CUBIC SPUN! 
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POTENTIODYNflMIC 50N05TLG 
Bc= 00.1579 Ba= 00.0399 




001 X SilOA 



67 



Figure 17. Tafel lines overlaid on PDP plot of Cu— Mn— Al— Fe— Ni alloy in 
synthetic seawater • 
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Figure 20. Cubic Spline derivatives of the anodic branch of Cu-Mn-Al 
alloy in synthetic seawater. 
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POTENTIODYNflMIC DLRLCLG 
TflFEL SLOPES; CUBIC SPUN 
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POTENTIODYNflMlC DL.RLCLG 
Bc= 00.0785 Ba= 00.0804 
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Figure 23. Tafel lines overlaid on PDF plot of Cu-Zn-Al 
synthetic seawater. 
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Figure 24. Cubic Spline derivatives of 630 series Bronze in synthetic 
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Figure 25. Tafel lines overlaid on PDP plot of 630 series Bronze in 
synthetic seawater - 
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POTENT I ODYNflMIC 017075 
Bc= 00.0383 Ba= 00.0116 
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Figure 27. Tafel lines overlaid on PDP plot of 7075 Aluminum alloy in 
3.5*/. NaCl solution. 
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POTENTIODYNRMIC SST304 
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Figure 29. Tafel lines overlaid on PDP plot of 304 Stainless steel in 
3.57. NaCl solution. 
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POTENTIODYNflMIC VCMOPD 
Bc= 00.1053 Ba= 00.2891 
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Figure 31. Tafel lines overlaid on PDP plot of Fe— Cr Mo alloy in 3.5% 
NaCl solution • 
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Figure 32. Cubic Spline derivatives of Fe-Cr-Al alloy in 3-57. NaCl 
solution • 
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Figure 33. Tafel lines overlaid on PDF plot of Fe— Cr A1 alloy in 3.5X 
NaCl solution. 
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Figure 34. Cubic Spline derivatives of Cu-Mn— A1 alloy in 3.5Z NaCl 
solution . 
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Figure 35. PDP plot of Cu-Mn-Al alloy in 3.5% NaCl solution. 
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Figure 37. Taiel lines overlaid on PDP plot of Cu-Mn-Al al 
NaCl solution. 
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Figure 38. Cubic Spline derivatives of 630 series Bronze in 3.5% NaCl 
solution . 
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APPENDIX B 

PAR 351 TO Z-248 DATA TRANSFER PROGRAM 



10 J = O 

20 SCREEN O 

30 COLOR 4,1,1 

40 KEy OFF:CLS:CLDSE 

50 DEFINT A-Z 

60 LOCATE 25,1 

70 PRINT STRINB$(60," ") 

80 FALSE = 0:7RLE=NDT FALSE 
90 X0FF?M<>^(19):X0N*=CI-R$(17) 

100 LIXT^TE 25,20:PRINT "PAR 351 TO Z-248 DATA TRANSFER PROGRAM" 

110 LOCATE 3,20:PRINT "YOU MUST PROPERLY CO^FIGLR TVt PAR 351." 

120 LOCATE 5,20:PRINT "AFTER SETTING THE TIME AMT DATE ON TVE 351" 

130 LOCATE 7,20:FRINT "PRESS 'SYSTEM MANAGER'" 

140 LOCATE a,20:FRlNT "PRSS COtvFIGLR SYSTEM " 

150 LOCATE 9,20:FRINT "SET 'PROTTOCDL' TO PRINT" 

160 LOCATE 10,20:PRINT "SET 'PARITY' TO MJ4E" 

170 LOCATE 11,20;PRINT "SET 'BALD TOTE' TO 2400" 

180 LOCATE 12,20:PRINT "SET 'STCP BIT' TO 1" 

190 LOCATE 13,20:PRINT "SET 'UDRD LENGTH' TO 8" 

■200 COLOR 14,1,1 

210 LOCATE 15,20:PRINT "PRESS ANY KEY TO CONTINUE" 

220 GO* = I(M<EY*:IF GD*= "" GOTO 220 
230 CLS 

240 COLOR 4,1,1 

250 LOCATE 25,20;PRINT "PAR :S5L TO Z-248 DATA TTONSFER PROGRAM" 

260 LOCATE 3,20:PRINT "PRE^ 'MAIN TEMJ' ON TFE 351" 

270 LOCATE 4,20:PRINT "PRESS 'FECALL/DISPLAY EXPERIfENT' " 

280 LOCATE 5,20:PRINT"PRESS 'DISPLAY EXPERIfENT' CN TTE COPY FROM DISK' 
SECTION" 

290 LOCATE 6,20:PRINT "SELECT TRE EXPERIfENT TO BE TFJANSFERED" 

300 LOCATE 7,20:PRINT "PRESS ' PLDT/FORMAT ' " 

310 LOCATE a,20:PRINT "ENSURE TFE X AXIS IS IN TFE LIfEAR SCALE" 

320 COLOR 2,1,1 

330 LOCATE 9,20:PRINT "IF NOT IN A LlfEAR DISPLAY TT-EN PRESS" 

340 LOCATE 10,20:FRINT " FORMAT DISPLAY, LlfEAR' , TEVIEW DISPLAY' , PLOT 
FCRMAT • " 

350 COLOR 14,1,1 

360 LOCATE 14,20:FRINT "PRESS ANY KEY TO CONTINUE" 

370 GO!t=IM<EY$:IF G0:*= "" GOTO 370 
380 CLS 

390 COLOR 4,1,1 

400 LOCATE 25,20:PRINT "FAR 351 TO Z-248 DATA TTW&ER PROGRAM" 

410 SPEED* = "2400" 

420 COfFIL* ="CDM2:"+SPEED*^",N,8,1,CS,DS" 

430 OPEN COfFIL* fG #1 
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"PflR 351 TO Z-248 DATA TRANBFER FTOSWI* 
"THIS PF03RRM lEES A RAM DISK TO SFtED 



"DO NOT ENTER Th€ DISK DRIVE CR T>E FILE TYPE" 
"HJT ON_Y TVE NATE CF TFE OUTPUT FILE" 



■ENTER TPE OUTPUT FILE NATE" 



4^ PRINT ttl, XCJFFt; 

450 (PEN FOR OUTPUT AS »2 

460 LOZATE 10,25:PRINT "(R)eceive a file OR (E)xit " 

470 LDCAIE 11,37: INUT TXRX$ 

480 IF (TXRX*0"R") AM) (TXRX^O"E") TVEN 460 
490 IF TXRX$="E" TPEN 1350 
500 CLS 

510 aXJOR 4,1,1 

520 LOCATE 25,20:PRINT "PAR 351 TO Z-248 DATA TTWEFER PROGRAM 
530 COLOR 15,1,1 

540 LOCATE 13,25:PRINT "FPESS 'PRINT DATA' ON TVE 351 TIEN" 

550 LOCATE 15,30:PRINT "PFESS ANY KEY TO CONTIMJE" 

560 GO* = IM<EY»:IF G0*="" GOTO 560 
570 CIS 

500 COLOR 4,1,1 
590 LCXT^TE 25,20:PRINT 
6CO LOCATE 10,20:PRINT 
EXECUTION. 

610 COLOR 2,1,1 
620 LCX:ATE 13,20;PRINT 
630 LOCATE 14,20:PRINT 
640 COLOR 14,1,1 
650 LOCATE 16,20;PRINT 
660 LCX^TE 18,25: INPUT FIL^ 

670 DSK!t="D:":DGKFILS=DSK$4FIL» 

680 OPENT "D:TTP" FCR OUTTUT AS #3 
690 N = O 
700 CLS 

710 PRINT #1, XON*; 

720 IF L0C(1)=0 THEN GDSLB 820 

730 IF LDC(1)>12B TPEN PAUBE=TPajE: PR I NT #1,X0FF* 

740 A*=INPUr*(l,#l) 

750 IF (A*=0-R*(10) ) OR (A*=OF»(13)) TTEN B*=OE»(13) :J=J+1 
IF (A(K>d«*(13) ) AND (A!K>OF{*(10)) TPEN PRINT #3,A$; 

IF (J > 3) AND (B«=CPRt(13)) AND (A*=OP»(13)) TPEN PRINT 
#3,CPFWi(13) 

700 IF (J > 3) AND (B*=CPm(13)) AND (Att=CPR^(13) ) TPEN B*=" ":NT=NH 
790 IF L0C(1)>0 TPEN 730 

000 IF PAUSE TPEN PAUGE=FALSE:PRINT 41,XCPI^; 

810 GOTO 720 

020 FDR 1=1 TO 5000 

830 IF LOCCDOO TPEN I = 9999 

840 NEXT I 

850 IF I>9999 TPEN RETURN 
860 CLOSE #3:CLS 
870 F*="D:"+A* 

800 OPEN "D:TPP" FOR INPUT AS #4 
890 CPEN DGKFILt FCR OUTTUT AS #3 
900 FOR I = 1 TO N 
910 X* = INPUT^(44,#4) 

920 Y* = INPUT*(2,#4) 



760 

770 
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930 B$ = 

940 C$ = 

9f» D% = Mim(B^,2,l) 

960 E*^ = MID^(B^,3,1) 

970 =: «ID^(B^,4,1) 

900 = MID$(B*,5,i) 

990 m = MID^(3$,6,1) 
lOOO = MID^(B^,7,1) 

1010 J% = MII>^(B^,0,1) 

1020 = MID^(B^,9,1) 

1030 FX$ = **^.000000’‘ 

1040 IF (E% = ”E”) TFEN MII»( X^, 18,9) = FX^ 

1050 IF (F^ = ••E“) TV€N MID$(X^, 18,9) = FX^ 

1060 IF (G^'*E") TFEN MID^( X^, 18,9) = FX$ 

1070 IF (H^‘*E**) TFEN MID^( X^, 18,9) = FX^ 

1000 IF (H*^”E") TFEN MID«( X^, 18,9) = FX$ 

1090 IF d^’*E”) TFEN MID^(X^,18,9)=FX^t 
lion I - (J^*’E^‘) TFEN MID^(X^,18,9)=FX^ 
iliO IF (K^‘*E”) TFEN MID1i( X^,18,9)=FX^ 

1120 F2^ = ’*E-12 ” 

1130 = MIE>^(X1i,42,l) 

1140 A1‘S=MID^(X^,43,1) 

1150 A2^II>t(X^,34,l) 

1160 A3^ = A^+^l^ 

1170 IF (A3^ = ■•12”) (A2^ = '•-” ) TFEN MID^( X^,39,5) = F2$ 

1100 = «II»(X^,1,1) 

1190 IF (P^ = ”P”) TFEN PRINT#3,X^ 

1200 rEXT I 
1210 CLOSE 413:CLS 
1220 J = 0 
1230 COLOR 1,5,4 
1240 BEEP 

1250 LOCATE 12,30:PRINT ”«$fTW|seFER COrPLETE*<<*” 

1260 LOCATE 13,30:PRIMr ”PRE^ ANY KEY TO CCMFINLE** 

1270 GD$=1M<EY^:1F GO^”” GOTO 1270 
1280 PRINT #1, XCFF^ 

12"?0 CLOSE #1 
1300 CLOSE #2 
1310 CLOSE nz 
1320 CLOSE #4 
1330 COLCF 1,1,1 
1340 GOTO 300 
1350 EM) 
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APPENDIX C 
CORROSION PROGRAM 



THIS IS THE MAIN PROGRPM. WITHIN IT H-E USER ENTERS Tl-E FILE 
NAME WHICH l-E/ShE WISFES TO ANALYZE. FOLLOWING THIS A ^E^4J 
(^=PEARS ALLOWING TFE OEER TO SELECT TVE DESIFED OFTICN. 



TFE FOLLOWING VWIBLES ARE CEFI^ED 



□DATA : 

CDATAl: 
MDATA : 

XDAT : 

YDAT : 

XDATl : 

YDATl : 

XI : 

X2 : 



XY 

Y1 

NB 

B2 

B1 

M 

I 

N 

01 

Ml 



J 

□2 

J2 

SI 



TVE 2 DirENSIONAL AFRAY FDR TFE DATA FILE CO^rrAINI^E 
POTENTIAL AND CURFENT DENSITY 

AS ABOVE BLIT OM-Y USED MEN OVERLAYING TWO CLJRVES 
A TWO DirENSIONAL ARRAY OF CDATA BLTT NOW IN TERE OF 
POTENTIAL AM) LOGIO CLRRENT DENSITY 
A OTE DIIENBICTMAL ARRAY CCNTAINirve THE X VALLES 
OF THE FIRST CURVE TO BE GRAPFED OR PLOTTED 
A 0^E DirENSIONAL ARWkY CONTAINING TVE Y VALLES 
OF TVE FIRST CURVE TO EE GRAPFED OR R.OTTED 
A OrE DirEhEIONAL AFRAY CONTAINING TVE X VALLES 
OF TVE SECOM) CURVE TO BE GRAPFED OR PLOTTED 
A C^E DltENSIOTMAL AFRAY CONTAINING TVE Y VALUES 
OF TVE SECOND CURVE TO BE GRAPFED OR PLOTTED 
TVE SLM OF TVE LOG OFRENT DE^BITIEB. USED IN TVE 
LIFEAR REGRESSION FDR DETERMINirsB TVE TAFEL CONSTANTS 
TVE SUM OF TFE SQUARS OF TFE LOG CURRENT DENSITY 
INITIALLY USED IN TFE MANNER A XI BL/T OTEE TFE 

TAFEL CONSTANT FIAS BEEN DETERMirED USED AS TFE LOG 
CURFENT DE^eITY VALUE IN CREATING TVE TAFEL LIIES 
TFE SLM OF TFE FRDDUCTS OF LOG CLRENT DETCITY AM) 
POTENTIAL. USED IN TFE SAFE MANER AS XI AM) X2 
TFE SLM OF TVE POTENTIAL. USED AS XI,X2,XY 
TFE NUTBER OF POINTS USED TO QEVERATE TFE TAFEL SLOPE 
TFE TAFEL SLOPE AS CALCULATED IN TFE REGRESSIGN 
TVE INTEFEEPT OF TVE TAFEL LItNE 
MENU FESPONSE 
DO LOOP COUNTER 

TVE INUBER OF POItTIS IN TVE DATA FILE HEirE FEAtf) 

TVE fNUBER OF POIf^S IN TVE FIRST DATA FILE, UNCHAM3ING 
TVE NUTBER OF POINTS IN TFE FIFST DATA FILE. USED IN 
TFE SUBFEUTirvE SLDFEl AM) CAN BE FETLF8NED WI TH A 
DIFFERENT VALLE 

TVE LOCATION OF ECDRR WITHIN M)ATA 
SATB AS J BUT F>ASSED AS AN AFEUEVENT TO SLOPEl 
SATE AS 01 BLfT USED OM_Y MEN O^EFLAYirvC DATA FILES 
SATE AS J BUT USED OM_Y MEN CNEFEAYINB DATA FILES 
TFE POSITION WITHIN MDATA MERE TFE FEOUIFETENTS OF 
TFE CENTRAL DIFFERNCE fETFED FAIL. CATHODIC BFJANCH 
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t 

t 

« 

t 

t 

t 

t 

t 

t 

t 

t 

t 

t 

t 

t 

t 

« 

t 

t 

t 

t 

« 

t 

t 

« 

« 

$ 

t 

* 



52 : 

53 : 



S4 

Ai 

A2 

A3 

A4 

Cl 

C2 

C3 

C4 

N0 



NWE : 
NPTE2 : 
TITLES: 

TITLE3: 
D1 : 

D2 : 

D3 : 
D4 : 
LIST : 



SATE AS SI ajT FOR TFE PMDDIC 

TVE TOTAL NUhOER OF POINTS IN TFE CATHODIC BF^^vCH 
WHICH FEET THE REOUIREFENTS FOR LEINS TFE CENTRAL 
DIFFEFENCE fETFIDD 

SATE AS S3 BLnr FOR TFE ANODIC BRANCH 

TFE POSITION WITHIN XDATl AM} YDATi WHICH MARKS TFE 

BEGINNING OF TFE TAFEL FEGION 

SAFE AS Al BUT TFE EM) OF TFE REGION 

SERV^ TFE SAFE PURPOSE AS Ai BLfT WITH RESPECT TO 

FDATA 

SAE AS A3 BLfT AS A2 

SA«E AS Al BLfT WITHIN XDAT AM) YDAT (CATHODIC) 

SAFE AS Cl BLfT TFE EM) OF TFE REGION 

SERVES TFE SAFE PURPOSE AS Cl BUT WITH RESPECT TO 

FCATA 

SAE AS C3 BLfT AS C2 

SET TO O. PASSED TO GRAHi AND PLOT TO ANOID TFE USE 
OF TFE EFPTY XDATl AFT) YDATi ARRAYS WFEN TIE ARE MDT 
USED 

TFE NAE CF TFE MAIN DATA FILE IN USE 

TFE NAE OF TFE SECOND FILE WFEN OERLAYING 

TFE SECOND HEADING FOR TFE GRAH As© PLOT. CHAEE5 

DEPEMENT ON TFE OPTION SELLCiHD 

RESERVED FOR FUTURE USE IF DESIRED 

CHARACTER STRING USED TO FOFJi TITLE2 

TFE CHARACTER FEFFESENTATICN OF B2 WFEN B2 REPRESEN 

TFE CATFCDIC TA=EL SLOPE 

CHARACTER STRING USED TO FORM TITLE2 

SAFE AS D2 BLfT FOR TFE ANODIC TAFEL SLOPE 

CHARACTER STRING USED FDR FESSAGES AS© PROFFTS 



REALM XDAT(2000),YDAT(20CX)), XDATl (2000), 

+YDATK2000) 

REAL»8 F©ATA(2000,2) ,MDATA1 (2000,2) ,X1 ,X2,NS,XY, XI ,Y1 
+ ,CDATA(2000,2) ,CDATA1 (2000,2) 

INTEGBR»2 I , J ,N,A1 ,A2,C1 ,C2,01 ,NI , J1 ,S1 ,S2,S3,S4,02, J2,NB 
INTEGERf2 A3,A4,C3,C4 

OTARACTERtO NP#E,NAFE2, TITLE2«2,YNA^»5,LI5T*50,TTTLE3*50 
CHARACTERtB D2,D4,D1*5,D3*6 
Tin_E2 = ' ' 

YNA^ = ' ' 

C 

C TFE SCREEN IS CLEARED AM) TFE FILE NAFE IS ENTERED. TFE USER 

C CAN ALSO EXIT THE PROGRAM OR GEl A LIST OF EXISTING DATA FILES. 

C THIS FRDFFT APPEARS IN PURPLE WHICH IS TFE COLOR USED WPEFEVER 
C A CHARACTER RESPONSE IS EXPECTED TO A PROFFT. IF A FILE NAFE IS 

C ENTERED WHICH DCES NUT EXIST AN ERFOT FESSAGE WILL APPEAR IN RED. 

C 

10 CALL QSMDDE(3) 

CALL QCL£AR(0,5) 
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15 



Lisr = ' ' 

CPU. QCMDV(15,12) 

U«rTE(6,0O5) LIST 
cm_ GCJtIV(25,12) 

LIST = ‘ ENTER TTE FILE NfiTE TO BE LGED' 

WRITE(6,0O5) LIST 
CPLL GOCV(20,11) 

LIST = A BLPM< TO EXIT, CF? •*DIR‘* FOR A DIRHJTCRY' 

WRITE(i.,805) LIST 
CPU- QCM>/(36,10) 

FEA0(6,800) Npre 
CPU- G01CV(40,9) 

IF (NPf« .EQ. ' ' ) TT€N 

GOTO 900 
EM)IF 

IF (NPfE .EQ. DIR ) THEN 
CPU. QCL£PR(0,5) 

CPLL QCM3/(1,23) 

PALEE ENTER ’’DIR *. /W" OR A BLPN< TO CCNTIMJE' 

CPU. QUSCFL(23,6,1,1,79,0,5) 

GOTD 15 
EM)IF 

CPU- DATAIN(CDATA,N,NPrC,IER) 

IF (lER .EQ. 1) TVEN 
GOTO lO 
ENDIF 

CPLL DATP«(CDATA,rt)ATA,N,J) 

CPU- OECK(rt)ATA,N,J,Sl,S2,S3,S4) 

20 CPLL QSMDDE(3) 

CPU- QCLEPR(0,2) 

J1 = J 

N = N-1 
01 = N 
N1 = N 
C 

C TFE MAIN TENU IS DISPLAYED IN GREEN. TIE COLOR GFEEN IS USED 
C WHENEVER PN INTEGER IS EXPECTED AS THE RESPONSE TO A PFOET. 

C IF TT-E IHPUT IS NOT PN INTEGER OR A OHPRACTER PN ERROR MESSAGE 

C IS DISPLAYED IN RED. 

C 

21 CPLL G01CV(15,22) 

LIST = ' ' 

WRITE(6,805) LIST 
CPLL QCMDV(15,18) 

LIST = 1. PLOT ORIGINAL POTENTIODYNPMIC CLFVE' 

WRITE(6,805) LIST 
CPU. QCM3/(15,17) 

LIST = '2. GEHERATE TPFEL SLOPES USING CUBIC SPLI^E METHOD 

MRITE(6,805) LIST 
CPU- QCM0V(15,16) 

LIST = '3. 6E^EFV^TE TPFEL SLOPES USING CENTFiPL DIFFEREME lETHDD' 



95 



14^ITE(6,805) LIST 
Ci=LL QCMDV( 15,15) 

LIST = '4. LEE flMJn-ER DATA FILE' 

WRITE(6,a05) LIST 

LIST = '5. OVER-AY TWO POTENTIODYNWIC CURVES' 

CALL QCMV(15,14) 

URITE(6,805) LIST 
CALL QCMDV(15,13) 

LIST = '6. EXIT' 

URITE(6,805) LIST 
CALL GDO^(29,12) 

LIST = ENTER 1,2,3, 4, 5 GR 6' 

WRITE(6,805) LIST 
CALL QCTCV(39,10) 

FEAO(6,810,BRR = 22,ICBTAT = J4) M 
CALL QCK]V(39,9) 

22 IF ((J4 .NE. O) .GR. (M .CT. 6)) Tt-EN 
C:ALL QEraDE(3) 

CALL QO.EAR(0,4) 

LIST = ' ' 

CALL QC>CV(20,17) 

URITE(6,805) LIST 

LIST = YGU DID NOT ENTER AN INTEGER GF' 

CALL QCMDV(20,16) 

URITE(6,805) LIST 

LIST = VALUE 1,2,3, 4 CK 5.' 

CALL QCMDV(20,15) 

WRITE(6,805) LIST 

LIST = 'FTESS ANY KEY TO CGNTINLE' 

CALL QCMDV(20,13) 

URITE(6,805) LIST 
CALL QCMO^(39,12) 

CALL QINKEY(I3,I4) 

GOTO 20 
EM)IF 
C 

C GPTIGN 1 HAS BEEN SELECTED. THIS GENERATES A GRAPH AfvO GPTIGNAL 
C RjOT CF Tl-E GRIGINAL POTENTIGDYNAfllC CURVE. 

C 

IF (M .EQ. 1) TTEN 
LO 25 I = 1,01 

XDAT(I) = M)ATA(I,2) 

YDAT(I) = MDATA(I,1) 

25 CCNTINJE 

N1 = N 
NB = O 

YNA»«E = VOLTS' 

TITLE2 = ' ' 

TITLE3=' ' 

CALL GRAPHl(XDAT,YDAT,01,Jl,NArE,XDATl,YDATl,NB,YNArE,TITLE2 
+ ,A1,A2,C1,C2,TITLE3) 
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eOTD 20 
EM>IF 
C 

C OPTION 2 BEEN SELECTED. THIS 6E^ERA7E5 7LE CERIVATIVEB Oh 
C THE POTENT ICDYNPmC CURVE LSINB TVt CUBIC SPLItvE MTIT-OD. 

C 

IF (M .EQ. 2) TTEN 
A1 = O 
A2 = O 
Cl = O 
C2 = O 
30 J1 = J 

Ml = N 

CPLL SL0FEl(r®ATA,Nl,Jl,NflME,Al,A2,Cl,C2) 

C 

C STILL WITHIN OPTION 2. TVE FRGGFW1 CHEEKS TO E^BLJRE THAT THE 

C TPFEL REGIONS FOR BOTH TVE ANODIC Afvtf) CATTCDIC FEGIGNB HAVE 

C BEEN SELECTED. 

C 

IF ((Cl .GE. C2) .OR. (A1 .GE. A2) ) TTEN 
CALL QSM]DE(3) 

CALL QCUEAR(0,4) 

LIST = ' ' 

CALL acn:N'(15,20) 

WRITE(6,0O5) LIST 

LIST = YOU DID NOT SELECT BOTH AN ANCDIC' 

CALL QOOV(15,19) 
l«RITE(6,805) LIST 

LIST = AM) CATHODIC BFVWCH OR ELSE TTE POINTS' 

CALL GCMDV(15,18) 

WRITE(6,0O5) LIST 

LIST = 'WERE NOT PROPEH_Y -mAFPED. IN ORDER TO' 

CALL GCM)V(15,17) 

WRITE(6,805) LIST 

LIST = 'PREVENT TTE PROGFWM FROM CRASHING YOU' 

CALL a:MDV(15,16) 

WRITE(6,0O5) LIST 

LIST = 'MUST START AGAIN.' 

CALL G3CMDV(15,15) 

WRITE(6,0O5) LIST 

LIST = PRESS ANY KEY TO CONTIKLE' 

C:ALL (3CMDV(15,13) 

WRITE(6,a05) LIST 
CTLl. QtM<EY(I3,I4) 

GOTO 20 
ENDIF 
C 

C ALL REDUIRTENTS HAVE BEEN SATISFIED. LIKEAR REEFESSICN IS 
C PERFORMED ON TEE TAFEL FEBIONS AM) TFE RESULTS ARE GRAPHED 
C 

D1 = Bc= ' 



97 



= 'M3LT5' 

[XI 35 I = 1,01 

XmJil) = M1ATA( 1,2) 

YDAT(I) = rt>ATA(I,l) 

35 QJVfTIMJE 

X2 = O 
XI = O 
XY = O 

NI = ABS(C2-C1) +1 
YI = O 

C3 = J-(4-K:1) 

C4 = J-(4-KI2) 

DO I = C4,C3 

XI = XI 4* nDATA(I,2) 

X2 = X2 + (rCATA(I,2)«2) 

XY = XY 4- (riDATA(I,2)m}ATA(I,l)) 

YI = YI+ MDATA( 1,1) 
oo cowrifvtJE 

NG = RJQAT(NI) 

Bi = ((XY*XI)-(YItX2))/((XI*XI)-(NB*X2)) 

B2 = ((XY*NS)-(YKXI))/((X2<N5)-(XI*XI)) 
Yl = M)ATA(J,1) 

XI = ((Y1-BD/B2) 

X2 = M)ATA(C4,2) 

DX = (X2-Xl)/(50.) 

DO 45 I = 1,50 

YDATl(I) = (B2<X2)4-B1 
XDATl(I) = X2 
X2 = X2- DX 
45 COrviTIMJE 

B2 = -B2 
CALL C0W(D2,B2) 

X2 = O 

XI = O 

XY = O 

NB = O 

YI = O 

A3 = J+4+A1 

A4 = J4-44A2 

DO 50 I = A3,A4 

XI = XI + m''TA(I,2) 

X2 = X2 4- (r«ATA(I,2)**2) 

XY = XY 4- (M)ATA(I,2)*M[)ATA(I,1) ) 

YI = YI4- MDATA(I,1) 

50 COViriNUE 

NI = ABS(A3~A4)4-1 
NB = FLQAT(NI) 

Bl = ((XY<XI)-(YI*X2))/((XI*XI)-(NB»X2)) 
B2 = ((XY<NB)-(YI*XI))/((X2*N5)-(XI*XD) 
Yl = MDATA(J,1) 
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XI = ((Y1-BD/B2) 

X2 = r®ATA(A4,2) 

DX = (X2-Xl)/(50.) 

DO 55 I = 51,100 

YDATl(I) = (B2*X1)+B1 
XDATl(I) = XI 
XI = X1+ DX 

55 CCNTIMJE 

D3 = ' Ba= ' 

CALL OXW(D4,B2) 

X2 = (XEATl(500)+XDATl(501))/2 
TITLE2=D1//D2//D3//D4 

^e = 100 

CALL GRP)PH1(XDAT,YDAT,□1,J1,NA^E:,XDAT1,YDAT1,NB,Y^4A^E,TI•^_E2 
+ ,Al,A2,Cl,C2,TnLE3) 

C 

C THIS PAFrr ALLOUe Tl-E LEER TO SELECT EITVER AMJD-ER AMIDIC CR 
C CATHODIC TAFEL REGION OR TO REVIEW TVE GRAPH WITHOUT HAVING 
C TO DETERMI^E BOTH TVE TAFEL REGIONS AGAIN. 

C 

56 CALL 0CLEAR(O,2) 

LIST = ' ' 

CALL GOl]V(20,16) 
l<«ITE(6,0O5) LIST 

LIST = 'IF TEE GRAPH DID NOT SHOW A GOOD CORRELATION' 

CALL 0CMDV(20,15) 

WRITE(6,0O5) LIST 

LIST = 'YOU CAN REPEAT THE CALCULATIONS FCR TEE' 

CALL OCMOV(20,14) 

WRITE(6,0O5) LIST 

LIST = 'BRANCH IN QUESTION' 

CALL GOfDV(20,13) 

WRITE(6,0O5) LIST 

LIST =1. REPEAT C«_CULATI0N5' 

CALL QCM]V(20,11) 
ltfRITE(6,805) LIST 
LIST = '2. EXIT THIS PORTION' 

CALL GChOV(20,9) 

WRITE(6,805) LISF 

LIST = ENTER A 1 OR A 2' 

CALL GOO/(20,7) 

C»RITE(6,805) LIST 
CALL QCM0V(25,6) 

READ(6,810,ERR=56) M 
IF (M .ED. 1) TEEN 
GOTO 30 
EM)IF 

IF (t1 .EQ. 2) TEEN 
GOTO 20 
El SF 

GOTO 56 
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c 

c 



EM) IF 
EM)IF 

TVE CENTWL DIFFEFNCE rtTHDD IS TG BE USED. TVE CFTICN PROCEEDS 
IN TFE SPTE MPNER flS CFTICN 2 

IF (M .EG. 3) TVEN 
A1 = O 
A2 = O 
Cl = O 
C2 = O 

57 CALL SL0PE(MDATA,Ol,Jl,N(W,Al,A2,Cl,C2,Sl,S2,S3,S4) 

IF ((Cl .GE. C2) .OR. (A1 .GE. A2)) TEEN 
(3SMODE(3) 

CALL QCLEAR(0,4) 

LIST = ' ' 

CALL OCM0W(15,2O) 

WRITE(6,805) LIST 

LIST = 'YOU DID NOT SOJEET BOTH AN ANCBIC' 

C7LL QCMDV(15,19) 

WRITE(6,805) LIST 

LIST = AND CATHODIC BRANCH OR ELSE TEE POINTS' 

CALL QCMDV(15,ia) 

WRITE(6,a05) LIST 

LIST = 'WEFE NJT mOPEFLY TRiiVEED. IN ORDER TO' 

CALL GCM:V(15,17) 

WRI1E(6,805) LIST 

LIST = 'PFEVENT TEE FTEOVM FROM OVEHING YOU' 

CALL QCM0W(15,16) 

U«ITE(6,a05) LIST 

LIST = 'MUST START AGAIN.' 

CALL QCM0W(15,15) 

WRITE(6,805) LIST 

LIST = PRESS ANY KEY TO CTIVITINLE' 

CALL GCMCM(15,13) 

WRITE(6,905) LIST 
CALL QIN<EY(I3,I4) 

GOTO 20 
EM)IF 
Dl='Bc= ' 

YNA^E = 'VOLTS' 

DO 60 I = 1,01 

XDAT(I) = W)ATA(I,2) 

YDAT(I) = M)ATA(I,1) 

60 CONTINLE 

X2 = O 
XI = O 
XY = O 

NI = ABS(C2-C1) +1 
YI = O 

C3 = J-(2+Cl) 
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C4 = J-(2+C2) 

CO 70 [ = C4,C3 

XI = XI + r«ATA(I,2) 

X2 = X2 + (M0ATA(I,2)**2) 

XY = XY + (rt)ATA(I,2)*«DATA(I,i)) 

YI = YI+ M3ATA(I,1) 

70 COMTIMJE 

re = FLOAT(NI) 

B1 = ((XY»XI)-(YI«X2))/((XI*XI)-(N5»X2)) 
H2 = ((XY«NS)-(YI*XI))/((X2»NB)-(XI»XI)) 
Yl = rOATA(J,l) 

XI = ((Y1-BD/B2) 

X2 = rCATA(C4,2) 

DX = (X2-Xl)/(50.) 

DO 00 I = 1,50 

YDATl(I) = (B2«X2)+B1 
XDATl(I) = X2 
X2 = X2- DX 
SO CCNTINLE 

B2 = -02 
CALL C0rV(D2,B2) 

D3 =' Ba= • 

X2 = O 

XI = O 

XY = O 

NS = O 

YI = O 

A3 = J+2+A1 

A4 = J+2+A2 

DO 90 I = A3,A4 

XI = XI + rt)ATA(I,2) 

X2 = X2 + (M>ATA(I,2)**2) 

XY = XY + (M)ATA(I,2)m)ATA(I,l)) 

YI = YI+ M3ATA(I,1) 

90 CONTINUE 

NI = ABS(A2-A1)+1 
NS = FLQAT(NI) 

B1 = ((XY»XI)-(YI»X2))/((XI»XI)-(NS*X2)) 
B2 = ((XY«ve)-(Yi*xi))/((X2*re)-(xi»xi)) 
Yl = rOATA(J,l) 

XI = ((Y1-BD/B2) 

X2 = MDATA(A4,2) 

DX = (X2-Xl)/(50. ) 

DO lOO I = 51,100 

YDATl(I) = (B2*X1)+B1 
XDATl(I) = XI 
XI = X1+ DX 
100 CONTINUE 

CALL C0rW(D4,B2) 

TIT1JE2=D1//D2//D3//D4 
NB = 100 
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58 



CflLL GRfimi ( XDAT , YDAT ,01 , J1 , NfiTE , XDATl , YDATl , NB,VN^yE , TITLES 
I- ,A1,A2,C1,C2,TITL£3) 

Cm_ QCL1£AR(0,2) 

LIST = ' ' 

C«-l- QD1CW(20,16) 

WRITE(6,805) LIST 

LIST = ' IF THE DID NOT 3-C3W A (TXT) CORRELATION- 

CALL QCMDN»(20,15) 

WRITE(6,0O5) LIST 

LIST = 'YOU CAW REPEAT Tl-E CALCULATIONS FOR TFE' 

CALL QCM0V(20,14) 

WRITE(6,805) LIST 
LIST = BRANCH IN GLE5TI0N' 

CALL QCMDV(20,13) 

WRITE(6,0O5) LIST 
LIST = '1. REPEAT CALCULATIONS' 

CALL QCMDV(20,11) 

V«ITE(6,0O5) LIST 
LIST = '2. EXIT THIS POROION' 

CALL QCMDV(20,9) 

WRITE(6,805) LIST 
LIST = ENTER A 1 OR A 2' 

CALL 0009(20,7) 

WRITE(6,805) LIST 
CALL (3009(25,6) 

READ(6,810,ERR=58) M 
IF (M .EQ. 1) TVEN 
GOTO 57 
EM)IF 

IF (M .ED. 2) THEN 
GOTO 20 
ELSE 

G(JTD 58 
EM)IF 
EM)IF 

C OFTICN 4 RETURre THE USER TO TFE BEGINMING CF TTE PROGRAM 

IF (M .EO. 4) TVEN 
GOTO 10 
EM3IF 

C CPTION 5 IS THE CNEFEAY. CN_Y TFE POrENTIODYNAMIC CLRVES CAW 
C BE GRAR-ED (XT PLOTTED. TFE (XTIGINAL FILE IS STILL A9AILABLE 
C F(JT FUFTTFER M^YMIRULATION 

IF (M .EQ. 5) 7FEN 
TirLE3=' ' 

105 CALL QCLEAR(0,5) 

LIST = ' ' 

CALL 0009(15,20) 



102 



750 



760 



800 

805 



U«ITE(6,805) LIST 
CPLL QCJCiO(25,12) 

LIST = ENTER ThC FILE TO BE USED' 

WRITE(6,805) LIST 
CALL 0Cm/(20,ll) 

LIST = 'A BLPM< TO EXIT, OR “DIR" FOR A DIRECTORY' 
URITE(6,a05) LIST 
CflLL QOCV(36,10) 

FEAD{6,800) NPTE2 
CALL QCMDV(40,9) 

IF (NATE2 .BQ, ' ' ) THEN 

GOTD 21 
ENDIF 

IF (NRI>E2 -ED. DIR') TT-EN 
CALL QCLEAR(0,5) 

CALL GOCW(l,23) 

P«JBE ENTER "DIR *. /W" OR A TO CONTINJE' 

CALL aUSCRL(23,6, 1,1,79,0,5) 

GOTO 105 
EM)IF 

DO 750 I = 1,01 

XDAT(I) = W3ATA(I,2) 

YDAT(I) = M>ATA(I,1) 

CONTIMJE 

CALL DATAIN(CDATAl,02,NAhE2,IER) 

IF (lER -ED. I) TVEN 
GOTO 105 
ENDIF 

CALL DATAM(CDATA1,MDATA1,Q2,J2) 

CALL CJ-ECK(M)ATAl,a2,J2,Sl,S2,S3,S4) 

DO 760 1=1,02 

XDATl(I) = r*t)ATAl(I,2) 

YDATl(I) = W)ATA1(I,1) 

CONTIMJE 

TITLE2 ='(J.€R_AY '//NAME//' ' //NATE2 
Y?^AfE = 'VOLTS' 

NA^E2=NA^E 
NATt = ' ' 

CALL 0RAPHl(XDAT,YDAT,01,Jl,NAft,XDATl,YDATl,02,YNArt,TITL£2, 
+ A1,A2, Cl, C2, TITLES) 

NATE = NATt2 
CALL aCLEAR(0,2) 

GOTO 21 
EM)IF 

IF (M .EQ. 6) THEN 
GOTO 900 
ELSE 

GOTO 20 
EM5IF 
FCRIAT(AB) 

F0R1AT(A5O,\) 
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810 FDRMAT(Il) 

900 CPLL GEMDDE(3) 

EM5 
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t THIS aiEFOJTIfvC; READS THE DATA FILE GEM3TTATED BY TFE 
» MODEL 351 WHICH HAS BEEN TFWNEFERFED CTHO A DISK 



C 

C 

C 

C 

C 

C 

C 



CDATA = TVE ARRAY WHICH CONTAINING POTENTIA. AN© 
ajRFENT DENSITY 

NATE = TFE NA-E OF TFE DATA FILE 

01 = CHARACTER STRING PRESENT IN TFE DATA FILE 

C2 = CHAfVCTER SrRI^G PRESENT IN TFE DATA FILE 

N = NCMBER CF DATA POINTS IN TFE FILE, 

DETEPm^ED DORire TFE PROGRAM 



SL®RDLfTirE DATAIN ( CDATA, N, NATE, lER) 

FEAL*8 CDATA(2000,2) 

INTEGER*2, N,I 

CHARACTER Cl»17,C2*6,NArE;»8,LIST*50,C3*50 

OPEN (UNIT = 10, ERR = 200, FILE = NATE, STATUS = OLD' ) 

DO 50 N=l,2000 

READ (10,900,ERR=60) C1,CDATA(N,1 ) ,C2,CDATA(N,2) 

50 CDNTINJE 

60 REWIND! 10) 

DO 70 I = 1,2000 

READ! 10, », ERR = 100) 

70 CONTIMJE 

100 IF !I .^E. N) TFEN 

CALL QSM0DE!3) 

CALL QCLEAR!0,4) 

LIST = ' ' 

CALL QCMCW!15,20) 

WRITE!6,905) LIST 

LIST = AN ERROR HAS BEEN ENCOLNTERED IN' 

CALL QCMDV!15,19) 

WRITE!6,905) LIST 

LIST = TFE FORMAT OF TFE DATA FILE. ' 

CALL QCM0V!15,18) 

WRITE!6,905) LIST 

LIST = THIS USUALLY INVOLVES A SWITCH BETWEEN' 

CALL QCMDV!15,17) 

WRITE!6,905) LIST 

LIST = 'FLOATING POINT AND EXPCTENTIAL NOTATION' 

CALL QCMDV! 15,16) 

WTITE!6,905) LIST 
LIST = ARDURD ECOOR. ' 

CALL QCMDV!15,15) 

U«ITE!6,905) LIST 

LIST = 'TO CORRECT THIS FHOBEM EXIT TFE PROGRAM' 
CALL QCMDV!15,13) 

WRITE!6,905) LIST 

LIST = ATvO USE TFE DOS ‘’EDLIN" FEATlfE. TFE' 
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D5LL GCMCiW(i5,12) 

WRITE(6,905) LIST 

LIST = PROBLEM BEGINS AT LI^e ' 

CALL 0CMUV(15,11) 
l««nE(6,905) LIST 
CALL QCMW(38,11) 
l«RITE(6,915) N 

LIST = 'PRESS ANY KEY TO CONTIMJE' 

CAUL QCM0V(15,8) 

WRI7E(6,905) LIST 
CALL QIM<EY(I3,I4) 

GOTO 990 
E^4)IF 
GOTO 990 

200 CALL C3SMDDE(3) 

CALL QCLEAR(0,4) 

LIST = • ' 

CALL QCMDL'nr,20) 

WRITE(6,905) LIST 

LIST = THE FILE NAI*E ENTERED DOES NOT EXIST.' 

CALL OCMOVI 15,19) 

WRITE(6,905) LIST 

LIST = TO EXAMINE TVE DIRECTORY OF POSSIBLE FILES' 
CALL QCM00(15,18) 

WRITE(6,905) LIST 

LIST = ENTER "DIR" M-EN ASKED TO ENTER TRE FILE' 
CALL QCM0Y(15,17) 

WRITE(6,905) LIST 

LIST = 'NA^E FOLLOWED BY THE DOS COTtlAND. ' 

CALL GCMD0(6,905) 

WRITE(6,905) LIST 
LIST = '"DIR *. /W"' 

CALL QCMDV(35,16) 

WRITE(6,905) LIST 

LIST = PRESS ANY KEY TO CONTIMJE' 

CALL QCMDV(15,14) 

U«ITE (6,905) LIST 
lER = 1 

CALL QIM<EY(I3,I4) 

GOTO 990 

900 FORMAT (A17,G8.4,A7,G10.4) 

905 F0RMAT(A5O,\) 

915 FORMAT! 14, \) 

990 RETURN 
END 
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t THIS SLER□Lm^E MPNIPULA7ES THE DATA IN TVE PEmv CDATA 
* THE PRCFER FTFM AND IDTENIFIES ECORR 



C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 



c 



c 



100 



CDATA = THE AmAY CONTAINING POTENTAIL PtO CURRENT 
DENSITY 

fDATA = CDATA COtsNERTED TO POTENTIAL AND LOG 
CURRENT DENSITY 

N = MJ1BER OF DATA POINTS IN T>E DATA FILE 

J = LOCATION WITHIN r®ATA WHERE ECORR IS LOCATED 

I = DO LOOP COLNTER 

XI = ABSOLUTE VALLE OF CURRENT DENSITY AT I 
USED FOR DETERMINING ECtJ« 

X2 = ABSOLUTE VALUE OF CURRENT DENSITY AT I-l 
USED FOR DETERMINING ECORR 
X3 = LGGiO(Xl) 

SUEFOJTINE DATAM (CDATA,rt>ATA,N, J ) 

REAL*8 MDATA(2000,2),CDATA(2000,2) 

INTEGBR*2 N,J 
X2 = 100000. 



COvNEFTT CURRENT DENSITY INTO LOG CURRENT DENSITY 

DO 100 I=l,N-2 

XI = ABS(CDATA( I ,2) ) 

X3 = LOGIO(XI) 

IF/TTEN TEETT UEED IN DETER1INING LOCATION OF ECOOR 

IF ( XI .LT. X2) TTEN 
J = I 
X2 = XI 
FT 

X2 = X2 
ENDIF 

MDATAd.l) = CDATA(I,1) 

MDATA(I,2) = X3 
CONTINLE 
RETURN 
EID 
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t THIS SLCROLfTIf^ PEFFCR-G TVD RJNCTIC3S6. ITS FIRST FUNCTIGN 
< IS TO CHECK TO ENSLre TT¥^T TVE PCJTEIsn-IPL VALUES FECORDED 
$ ACCURATELY REFLECT THE ACTUAL POLARIZATICM. TTE SECOND IS 
t TO DETERMIfvE TVE NUTBER OF POI^fTS IN THE AfsEDIC AM) CATHTOIC 

* BRAMKES WHICH MAY HE USED IN TVE CENTRAL DIFFERNCE rETHDD. 

« 

♦ MDATA 

« 

* N3 

♦ 

* N4 

« 

* OTHER VARIBLES ARE USED fEFELY IN OECKINB REQUIRED CONDITIONS 

« 

SU8FaJTI^E CFECK(r^TA,N,J,Nl,N2,N3,N4) 

REAL*8 M)ATA(2000,2) 

INTEGER<2 I ,J ,N,N1 ,N2,N3,N4 
CHARACnER*50 LIST 

C TFE FIRST INCRE^€NT IN POTENTIAL FDR TTE CATHODIC BRAMK IS 

C DETERMINED 

10 Y = MDATA(J,1) 

Y1 = MDATA(J-1,1) 

Z1 = ABS(ABS(Y1)-ABS(Y)) 

Z1 = ANINr(Zl/0.0001) 

C SOBSEQLENT POTENTIAL INCREMENTS, WHICH MiST BE EOL^L TO DE 

C FIRST ARE OECKED. LOG CURRENT DENSITY VmJES ARE ALSO OECKED 

C TO ENSURE THAT TEEY WILL NOT RESULT IN TEE CALCULATION OF A 

C SLOPE OF INFINITE VALLE. 

DO 20 I = 

Y = MDATA(I,1) 

Y1 = MDATA( 1-1,1) 

X = MDATA(I,2) 

XI = MDATA( 1-1,2) 

Z = ABS(ABS(Y1)-ABS(Y) ) 

Z = ANINT(Z/0.0001) 

C OTE OF TVE TWO REQUIRED COMDITIONS HAS FAILED. 

IF ((Z .^E. Zl) .OR. ( XI .LE. X )) GOTO 25 

20 CONTINUE 
25 N1 = 1+2 



: THAT ARRAY CONTAINING POTENTIAL AM) LOG CLRRENT 
DENSITY 

: TVE NUMBER OF POINTS IN TVE CATHODIC ERANCH WHICH 
CAN BE USED IN TVE CENTRAL DIVFEhtNCE VETHDD 
: TVE NUMBER OF POINTS IN TVE ANCBIC BRANCH WHICH 
CAN BE USED IN TVE CENTT^ DIFFERVJCE rETHDD 



lOG 



n n 



ThE FIRST INCRB-Ofl" IN POTtNTI^ FCR T>E ANCDIC BRfiNCH IS 
DETERIillVED 

X = M)ATA(J+l,i) 

Y = NE)ATA(J+2,1) 

Z1 = ABS(X-Y) 

Z1 = ANIMT(Z1/0.CXX)1) 

C SUBSeUEMT POTTENTIAL INCFOENFS, UHICH rUST BE EQLK^ TO THE 

C FIRST fiFE CHECKED. LOG CURRENT DENSITY VmJES ARE ALSO CHECKED 

C TO ENSURE THAT THEY WILL NOT RESULT IN THE CALCULATION OF A 
C SLCFE OF IIVFINITE VALUE. 

DO 30 I = J+2,N-1 

Y = r-CATA( 1,1) 

Y1 = M)ATA( 1+1,1) 

X = MDATA(I,2) 

XI = r®ATA( 1+1,2) 

Z = ABS(Yl-Y) 

Z = ANINT(Z/0.0001) 

C GTE OF THE TWO FEOUIRED COM)ITIONS HAS FAILED. 

IF ((Z .NE. Zl) .OR. ( XI .LE. X )) GOTO 35 

30 CONTINUE 

35 N2 = I - 2 
N3 = J - N1 
N4 = N2 - J 

C THE NUMBER OF HDINTS AWAILAELE FOR EACH BRANCH IS COEARED 
C WITH A MINUMIUM NUMBER OF 15. IF HEITHER BRANCH HAS ENOUGH 

C AVAILABLE POINTS VARIOUS HESSACES ARE DISPLAYED LISTING 

C POSSIBLE CAUSES AMD OPTIONS. 

IF ((N3 .LE. 15) .OR. (N4 .LE. 15)) THEN 
CALL QCL£AR(0,4) 

LIST = ' ' 

CALL QCMCN'(20,20) 

WRITE(6,805) LIST 

LIST = THE CENTRAL DIFFEF3EE fETHCD REQUIFES THAT’ 

CALL QCMDV(20,19) 

WRITE(6,805) LIST 

LIST = THE CURRENT BE MEASURED AT EQUAL DIFFERENTIALS' 

CALL QCMDV(20,1B) 

WRITE(6,805) LIST 

LIST = OF THE APPLIED VOLTAGE. THE PAR MXGL 351' 

CALL QCMDV(20,17) 

WRITE(6,805) LIST 

LIST = RECORDS THE CURRENT AT INCREHENTS OF 2.0 (IN'' 

CALL QCMDV(20,16) 
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400 



V«nE(6,e05) LIST 

LIST = 'UVITIL TVE IN CURRENT EXCEEDS A PRESCRIBED' 

COLL QDtJV(20,15) 

WRrTE(<b,a05) LIST 

LIST = LEVEL. AT TWT POINT Tl-E 351 BEGINS FECORDING' 
CALL QCm«'(20,14) 

WRITE(6,0O5) LIST 

LIST = 'TPE CURRENT AT INCREMENTS OF 0.5 ahf. AT AN' 

CALL GCMDV(20,13) 

WRI7E(6,805) LIST 

LIST = 'HIGPER APFEIED VOLTAGE TIE CHASEE IN CURRENT' 
CALL QCMDV(20,i2) 

WRITE(6,805) LIST 

LIST = 'EECREASES TO A POINT WHERE TIE CURRENT IS AGAIN' 
CALL QCMDV(20,11) 

WHnE(A,a05) LIST 

LIST = 'PEAGURED AT INCREPENTS CF TT-E APPLIED VOLTAGE' 
CALL QCMDV(20,10) 

WHITE(6,0O5) LIST 
LIST = EQUAL TO 2.0 mV' 

CALL QCM0V(20,9) 

WRnE(6,805) LIST 

LIST = 'PRESS ANT KEY TO OJNTINLJE' 

CALL QCPCV(20,7) 

WRITE(6,a05) LIST 
CALL QIW<EY(I,MB) 

CALL QCLEAR(0,4) 

LIST = ' ' 

CALL QCMDV(20,20) 

WRITE(6,a05) LIST 

LIST = TTHE DATA RECORDED FOR THIS EXPERIPENT El TVER : ' 
CALL QCraV(20,19) 

W^ITE(6,a05) LIST 
LIST = ' ' 

CALL QCM0V(6,18) 

WRITE(6,805) LIST 

LIST = 1. REFLECTS THAT TVE CHANGES IN CURRENT' 

CALL QCM0V(20,17) 

WRITE(6,e05) LIST 

LIST = MEAeURETENT DESCRIBED EARLIER DID NOT OCCUR' 

CALL GCMDV(20,16) 

WTITE(6,0O5) LIST 

LIST = AT A POINT(S), RELATIVE TO ECODR, THAT WOULD' 

CALL QCM0V(20,15) 

WRITE(A,0O5) LIST 

LIST = GENERATE ENOUGH DATA POINTS FOR TVE USE OF TVE' 
CALL QCM0V(20,14) 

UFITE(6,a05) LIST 

LIST = CENTRAL DIFFERENCE VETHOD. ' 

CALL QCM0V(20,13) 

UFITE(6,805) LIST 
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LIST = 'OR' 

CPLL 00*0/(39, 11) 

WRirE(6,0O5) LIST 

LIST = '2. TTE MODEL 351 DID NOT FtCOFlD TT-E APPLIED' 
CALL 000/(20,9) 

WRI7E(6,805) LIST 

LIST = POTENTIAL PBOPEFiLY. INIS CAN BE CCRRECTED' 

CALL 00*0/(20,0) 

WRITE(6,0O5) LIST 

LIST = • IF TVC PROBLEM IS NOT (XPRECTTED T>€N Tl-E’ 

C:ALL 00*0/(20,7) 

WRITE(6,805) LIST 

LIST = CENTRAL DIFFERNCE tETHCD WILL NOT WCF3<, AM) THE' 
CALL 0CMDV(20,6) 

WRITE(6,805) LIST 

LIST = 'CUBIC SPLI^E l«eTTHDD MAY NOT UCFK' 

00*0/(20,5) 

WRITE(6,0O5) LIST 

LIST = PRESS ANY KEY TO CCNTIMJE' 

CALL 00*0/(20,3) 

WRITE(6,0O5) LIST 
CALL QIM<EY(I,I*B) 

CALL (3CLEAR(0,4) 

LIST = ' IF AN ATTE1*PT IS MADE TO CCFFECT TIE' 

CALL 000/(20,17) 

U«nE(6,a05) LIST 

LIST = PRCBLEM TFE CCFRECT fETTHCD IS DETERMINED' 

CALL (300/(20,16) 

WRITE(6,a05) LIST 

LIST = ' IF TT-E ATTEJ*f=T IS OvEUCCESSFUL 71-EN TT-E' 

CALL 010*0/(20,15) 

WRITE(6,a05) LIST 

LIST = CUBIC SPLINE NETHDD SPDLLD BE USED.' 

CALL 00*0/(20,14) 

WRITE(6,805) LIST 

LIST = ' IF TPE ATTENPT IS SUCCESSFUL EITTER’ 

CALL (30*0/(20,13) 

WRITE(6,805) LIST 

LIST = 'NETH3D MAY BE USED.' 

C7LL 00*0/(20,12) 

UFITE(6,E305) LIST 

LIST = ' IF THIS SET CF NESSAGES APPEARS AIDIN' 

CALL (30*0/(20,10) 

WRITE(6,0O5) LIST 

LIST = AFTER ATTENPTING TO (XRFECT TPE PROBLEM' 

CALL (30*0/(20,9) 

WRITE(6,a05) LIST 

LIST = THEN TT-E (OJRRENT DENSITY VALLES ARE SUCH' 

CALL (3(00/(20,8) 

MRITE(6,0O5) LIST 

LIST = THAT NEITPER METHOD Cfti BE USED' 
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CPLL GOtDV(20,7) 

WRITE(6,805) LIST 

LIST = PRESS KEY TO CONTINUE' 

CPLL GOt3V(20,5) 

WRITE(6,805) LIST 
CPLL QIM<EY(I,MB) 

420 CPLL QCLEP«(0,2) 

LIST = '1. DG MDT CORRECT TTE PROBLEM.' 

CPLL QCM00(25,13) 

WRITE(6,805) LIST 

LIST = '2. IRY CORRECT TEE PROBLEM.' 

CALL QCMD0(25,11) 

WRITE(6,0O5) LIST 

LIST = ENTER A i OR A 2' 

CALL QCM0V(32,9) 

WRITE(6,0O5) LIST 
CALL QCM00(39,8) 

REA0(6,^T0,ERR = 450, lOSTAT = J4) M 
450 IF ((M .GT. 2) .OR. (J4 .fvE. O) ) TEEN 

CALL QCL£AR(0,4) 

LIST = ' ' 

CALL QCMD0(21,16) 

WRITE(6,0O5) LIST 

LIST = 'YOU DID NOT ENTER AN INTEGER OF' 

CALL QCM0O(2O,15) 

WRITE(6,B05) LIST 
LIST = 'VALUE 1 OR 2' 

CAUL GL2«V(37,14) 

WRITE(6,0O5) LIST 

LIST = 'PRESS ANY KEY TO CONTINUE' 

CALL QCMDV(25,12) 

URITE(6,0O5) LIST 
CALL QCMDV(39,li) 

CALL QirvKEY(I,MB) 

GOTO 420 
EEff)IF 

IF (M .ED. 1) TEEN 
N9 = 1 
GOTO 1000 
EM)IF 

C AN ATTEMPT IS BEIM3 MADE TO CORRECT TEE POTET'^AIL VALUES SUCH 
C THAT TT-EY ACCURATELY REFLECT THE EXPERirENTAL POTENTIALS. 

C ELOCKS 500 AND 600 CORRECT TIE PRCH-EM, BUT OTsLY ONE OF TFE 

C TVO ELOCKS ARE USED FOR ANY GIVEN FILE. 

IF (M .ED. 2) TFEN 
DO 500 I = 1,N,2 
Y1 = MDATA(I,1) 

Y2 = rDATA( 1^1,1) 

IF (Y1 .ED. Y2) TI-EN 
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500 



600 



805 

810 

1000 



l*®ATA(I,l) = MDATA(I,1) - 0.0005DO 
N9 = 1 
END IF 
CONTIMJE 
DO 600 I = 2,N,2 
Y1 = rCJATAd,!) 

Y2 = M)ATA( 1+1,1) 

IF (Y1 .03. Y2) THEN 
MDATA(I,1) =M>=kTA(I,l) - O.0OO5DO 
N9 = 1 
ENDIF 
CONTirUE 

IF (N9 .ED. 1) TT-EN 
CflLL QCLEfiR(0,2) 

LIST = ' ' 

Cm_ GCMD0(3O,14) 

WRITE(6,805) LIST 

LIST = COFFECTION SUCCESSFUL' 

CfiLL GCMD0(3O,13) 

WRITE(6,805) LIST 

LIST = 'USE EITHER l*ETHDD' 

CflLL QCMCIV(30,11) 

WRITE(6,805) LIST 

LIST = ' P RESS ANY KEY TO CONTINUE' 

CPLL GOCJ0(30,9) 

WRITE(6,B05) LIST 
CPtJL QIM<EY(I,MB) 

GOTO 10 
ELSE 

CfiLL QCLEAR(0,4) 

LIST = 'CORRECTION UNSUCCESSFUL' 

C«_L GCMDO(30,13) 

WTITE(6,805) LIST 

LIST = USE EITFER CUBIC SFLir>E hETHOD' 
CfiLL QCMa0(30,ll) 

WRITE(6,805) LIST 

LIST = 'PRESS ANY KEY TO CONTIMJE' 

CALL GCM0V(3O,9) 

WRITE(6,805) LIST 
CALL QIM<EY(I,MB) 

GOTO lOOO 
ENDIF 
ENDIF 
EM)IF 

FDR1AT(A50,\) 

FORMAT (II) 

END 
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t 

% 

* 

t 

t 

t 

t 

t 

t 



XDAT 

YDAT 

XDATl 

YDATl 

TX 



TY 

TXL 

TYl 



Tit X VPL.LE5 CF TfE FIRST GR^ 

Tl-E Y VPEUEB CF T>£ FIRST GR^H 
TTE X VftUES GF TTE SECIF® GRPPH 
THE Y WLLJES CF TFE SECTF® GRAPH 

XDAT AS CFIGINPE PASSED AfvD SCALED. USED TO ALLOW 

TTE U3EJR TO RESCALE TTt AXIS AT® mirTTAIN ALL TTE 

VALLE ORIGINALLY PASSED 

USED AS TX BUT FOR YDAT 

U5tD AS TX BUT FOR XDATl 

USED AS TX BUT FOR YDATl 



C IN THIS SUBROUTirt ALL PROFTS ARE FtFLIED TO BY TVE USE CF 
C SOFT KEYS. THIS MEA#t THAT A CARRIAGE RETUFFJ IS NOT USED. 



SUBFOJrirt GRAFH1(XDAT,YDAT,^^3,J,^¥^,XDAT.1 ,^^T1 ,N4,YNA^E,T1TLE2, 
1 , A2 , C 1 , C2 , T I TLE3 ) 

REAL»4 XDAT(2000) ,YDAT(20CO) ,XDAT1(2000) ,YDAT1(2000) ,F,F1,F2 
REALM TX (2000), TY(2000),TX1 (2000) ,TY1 (2000) ,X2( 1) , Y2( 1 ) 

IN1LGER<2 K2,K3,K5,K6,L,XL,XR,YL,YR,^e,Al,A2,Cl,C2,M,^A(80) 
IN1EGER*2 I,J,N,K7,K8,K9, JCOL(IO) , JROW(IO) ,N1,I®,N3,N4,N5,N6,N7 
IN1EGER*2 L1,L2,L3,L4 

CHARACTER YNArt*5,FACT$^,XNA<ttl2,TITLE*15,TITLE2*32 

CFmA£TERir28, QUE5, TITLE3*50 
N7 = 1 
M = O 

N = MAX(N3,N4) 

N3 = N3 - 3 

ISi4 = fs|4 - 3 

IF (N4 .LE. O) N4 = 1 

C TEE ARRAYS ARE SCALED SUCH THAT TVE Y AXIS VALUES ARE DiaiERNIBLE 
C FROM EACH OTIER 

F = riIN(ABS(YDAT(l)),ABS(YDAT(N3))) 

FI = MIN(^«S(YDAT1(1)),ABS(YDAT1(N4))) 

IF (N4 .GT. 2) TEEN 
F2 = MIN(F1,F) 

ELSE 

F2 = F 
ENDIF 

IF (F2 .GE. 1.) THEN 
F2 = 1.0 
FACT = 'x 1' 

EM)IF 

IF ((F2 .GE. .1) (F2 .LT. 1.)) TEEN 

F2 = 10 
FACT = X 10' 

ETOIF 
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IF ((F2 .GE. .01) .PM). (F2 .LT. .1)) TVGM 
F2 = 100. 

FPCT = X 100* 

EM)IF 

IF ((F2 .GE. .OOl) .PM). (F2 .LT. .01)) 71-EN 
F2 = 1000 
FPCT = X 1000' 

ELSE 

IF (F2 .LE. 0.0010) 71-EN 
F2 = 1000 
FPCT = 'x 1000' 

EM)IF 

ENDIF 

DO 6 I = 1,2000 

YDAT(I) = F2»YDPT(I) 

YDATl(I) =F2*YDAT1(I) 

6 CDNTIMJE 

5 N5 = M3 

N6 = N4 

C 7VE PEFJ1PIVCNT PFRAYS PRE FILLED 

DO 10 I = 1,2000 
TX(I) = XDPT(I) 

TY(I) = YDAT(I) 

TXl(I) = XDATl(I) 

TYl(l) = YDPTl(I) 

10 CONTIMJE 

C D-E GRPPH AXIS PND WIM)GW PFE DETERlirED 

XmX = -1000 
XMIN = 1000 

YMAX = xmx 

YMIN = XMIN 

XMAXl = XMAX 

XMINl = XMIN 

YMAXl = XMAX 

YMINl = XMIN 

JCOLl = 90 

JC0L2 = 610 

JROWl = 35 

JFCW2 = 300 

DO 200 I = 1,N3 

YMAX = MAX(YDAT<I),YMAX) 

XMAX = MAX(XDAT(I),XmX) 

XMIN = MIN(XDAT(I),XMIN) 

YMIN = MIN(YDAT(I),YMIN) 

200 CONTIMJE 

205 DO 210 I = 1,N4 

XMAXl = MAX(XDAT1(I),XMAX1) 
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YMAXl = mX(YDATl( I) ,YmXl) 

XMINi = liIN(XDATl( 1) ,XMIN1) 

VniNl = hIN(YDATi(I),YMINl) 

210 CXJvrriMJE 

215 YFIN = (MAX(YmX,YmXl) ) + .1 
YST = (MIN(YMIN,YMIN1) ) - -1 
XST = (MIN( XT1IN, XMINI ) ) -.1 
XFiN = (r^x(xmx,xMAXi) ) + .1 
225 YMIN = YST 
YMftX = YFIN 
XmX = XFIN 
XMIN = XST 
YORG = YMIN 
XCFG = XMIN 
XfvWTE = 'LOG I A/cm2' 

TITLE = PUYEmiam^miC' 

ICPT = O 
YO^ERX ^ 0.3 
ASPECT = 1.5 

C TFE GRAPH AXIS ARE DRAWM AM) LABELED. TITLES ARE PLACED 
CALL QSMjDE( 16 ) 

CALL GPLDT(JCC1JL,JCGLJ2,JR0W1,JRDW2, XMIN, XMAX, YMIN, YMAX,XCPG,YDRG, 
+ICFT,Yr]OERX , ASPECT) 

XmiOR = (XFIM-XST)/10. 

XMINCR = O 
LABEL = -1 
M)EC = 2 

YMAOCF =(YFIN-YST)/10. 

YMINDR = O 

rsKNT = 5 

MXJTS = O 

IKOLOR = 1 

ISY^B□L=-2 

KLRSYM=1 

CALL QC0LCF(O,O) 

CALL QSETlJP(MX^■S,II<^JDR,ISY^01,,KL^ 

CALL QXAXIS(XST, XFIN, XMAJOR, XMINCR, LABEL,M)EC) 

CALL aYAO<IS(YST,YFIN,YMAOCJR,YMI^IT?,LAflEL,^OEC) 

CALL QXTICS(M<Nr,JCCL,JRC]W) 

CA4J^ GGTXT(12,XNAriE,2,340,l,0) 

CALL aYTICS(NKNr,JCOL,JRGW) 

CALL QGTXT ( 5, YNAME, 2, 1,246,-1) 

CALL GGTXT(6, FACT, 2,1,162,“!) 

CALL QGTXT( 15, TITLE, 2, 254, 335,0) 

CALL GGTXT(8,NAI^, 2,382, 335,0) 

CALL QGTXT(32,TITLE2,2,254,321,0) 

CALL QGTXT(50,TITLE3,2,182,307,0) 

IKOLOR = 2 

CALL QSETIJP(MXJT5,IK0LCJR, I SYMBOL, KLRSYM) 
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CPU. C3RTDI(XST,YST,L1,L2) 

CPLL CFTOI(XFIN,YFIN,L5,L4) 

CPLL 0Lrive(Ll,L2,Ll,L4,2) 

CPLL (]LINE(L1,L4,L3,L4,2) 

CALL a_I^E(L3,L4,L3,L2,2) 

CPU. 0_INE(L3,L2,L1,L2,2) 

C TVe FIRST CURVE IS DRPWN. Tl-E LINE CCLCF IS 0«3N3ED m) TVE 
C SECOM3 CLFVE IS DRAWN 

CALL OTAa_(l,N5,XDAT,YDAT) 

IKd-CR = 1 

CALL QSETlJP(^mTS,II<a_(J^,ISY^«CL,KLRSY^I) 

CALL GTAlBL(l,Nfe,XDATl,YDATl) 

CALL QBEEP 
M = M + 1 

C IF TVE CURVES ARE OF DERIVATIVES A LirEAR REGICN HAS ALJEADY 
C BEEN SELECTED, TVE USER MAY DECIDE TO GE^ERATE Tl-E TAREL CONSTANTS 
C BAEED ON THIS PORTION OF THE CURVE. 

300 IF (( M .GT. 1) .AND. (N4 .GE.O) ) TVEN 

QUE5 = 'GEhEFWTE TAFEL CONSTANT- 
CALL QGTXT(28,GLES,3, 1,335,0) 

GOES = ' Y OR N ?' 

CALL GGTXT(20,OUES,3, 1,320,0) 

CALL GB F EP 
CALL QIN<EY(I,M8) 

IF ((MB .EQ. 89) .OR. (MB .EQ. 121)) TVEN 
GLEB = ' ' 

CPLL QGTXT(28,QUES,3,1, 335,0) 

CPU. GETXT(28,QUE5,3, 1,320,0) 

3B1 (3UES = ' IS THIS TVE ANCDIC OR' 

CPLL GSTXT(28,QUES,3, 1,335,0) 

OJES = 'CPTVCDIC BRANCH- 
CALL GUrXT(28,(aJES,3, 1,320,0) 

(IFS = - A CF C ?- 
CPLL GSTXT(28,QUES,3, 1,305,0) 

CPLL GEEEP 
CPLL QIM<EY(I,MB) 

IF ((MB .EQ. 65) .OR. (MB .ED.97) ) TVEN 
A1 = KB 
A2 = K9 
GCTTO 390 
EM)IF 

IF ((MB .EG. 67) .OR. (MB .EG. 99)) TVEN 
Cl = K6 
C2 = K7 
GOTO 390 
ELSE 

GOTO 381 
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c 

c 

390 



400 



410 

420 



EM}IF 

EM)IF 

IF ((r« .EG. 78) .OR. (MB .EG. 110)) TVtN 
GQrO 390 
ELSE 

GOTO 380 
ENDIF 
EM)IF 

BY TFH USE OF TVE OJF2SOR KEYS ON TEE MJTERIC KEYPAD TEE USER 
SELECTS A NEW SET OF AXIS TO BE PLOTTED 

GEES = ' ' 

CALL aGTXT(28,0LES,3,l,335,0) 

OPEL Q8TXT(2B,QLE5, 3, 1,320,0) 

CALL □GTXT(28,QLE5, 3, 1,305,0) 

QLES = r^EW AXIS Y OR N ?' 

CALL GGTXT(28, QLES, 3, 1,335,0) 

CALL GBEEP 
CALL QIN<EY(I,MB) 

IF ((MB .EG. 89) .OR. (MB .EG. 121)) TEEN 
GOTO 400 
ELSE 

IF ((MB .EG. 78) .OR. (MB .EG. 110)) TEEN 
GOTO BOO 
ELSE 

GOTO 390 
ENDIF 
EEOIF 
GOES = ' ' 

CALL GGTXr(28, GOES, 1,1, 335,0) 

C:ALL GGTXT(28,GUES, 1,1,320,0) 

OLES = '^EW,LAST,(^T (3TIGINAL (N,L,D)' 

CPLl. GGTXT( 28, DUES, 3, 1,335,0) 

CIALL GBEEP 
CAIX QIN<EY(I,MB) 

IF ((MB .EG. 79) .OR. (MB .EG. Ill)) TEEN 
DO 410 I = 1,N3 
XDAT(I) = TX(I) 

YDAT(I) = TY(I) 

COrTTINLE 
DO 420 I = l,IM4 
XDATl(I) = TXl(I) 

YDATl(I) = TYl(I) 

CONTINJE 
M = O 
GOTO 5 
ENDIF 

IF ((MB .EG. 76) .OR. (MB .EG. 100)) TEEN 
XST = XSTl 
XFIN = XFINl 
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430 



500 



YST = YSTl 

= YFINl 
N5 = he 
N£» = M6 
GOTO 530 
EM>IF 

IF (( MB .EQ. 78) -CR. (MB .EQ. 110)) 71-EN 
XSTl = XST 
YSTl = YST 
XFINl = XFIN 
YFINl = YFIN 
M5 = N5 
M6 = N6 

(3UES = 'LEE T>E OJRSCR KEY TO M0S«E' 

CflLL GGTXT(28,GUES, 3, 1,335,0) 

QLES = TO n-E LOUER LEPT^ 

CfiLl_ (XrrXT(28,GLES,3, 1,320,0) 

CUES = CIFJvER, HJE: TO EXIT' 

CfUL QGTXT(28,ajES,3, 1,305,0) 

CALL GEEEP 
ClfiLL QLSTAT(l) 

CIALL GLTYPE(3,3) 

CALL □LMCN'(320,175) 

CALL QIN<EY(I,L) 

IF (L -ED. 72) n-EN 
C:ALJ_ (3LM)WR(0,5) 

(XJTO 500 
ENDIF 

IF (L .EQ. 75) TT-EN 
CALL GLM0WR(-5,O) 

GOTO 500 
ENDIF 

IF (L .EQ. 77) TVEN 
CyiLL OJ1J^(5,0) 

GOTO 500 
EM)IF 

IF (L .EQ. 80) TTEN 
CALL (3_hCMR(0,-5) 

GOTO 500 
ENDIF 

IF (L .EQ. 71) TVEN 
CALL QLFEAD(XL,YL) 

CIALL QLMCNR(320-XL,175-YL) 

CALL QITOR(XL,YL,XST,YST) 

X2(l) = XST 
Y2(l) = YST 
CALL QSETIJP(5,3,43,3) 

CALL QTAa_(0,N7,X2,Y2) 

GOTO 510 
ELSE 
GOTO 430 
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510 



520 



530 



EM)IF 

GLE5 = 'U5E TFE OJRSCF KEY TO 
CPLL QGTXT( 28, CUES, 3, 1,335,0) 

QUES = 'TO THE LFPER RIGHT 
CPU. QGTXTCS,QUES, 3, 1,320,0) 

QLES = 'CORSER, HDhE TO EXIT' 

CALL QGTXT(28,QUES, 3, 1,305,0) 

CALL WEEP 
CALL QIM<EY(I,L) 

IF (L .EQ. 72) THEN 
CALL GLJiCMT(0,5) 

GOTO 520 
EM)IF 

IF (L .EQ. 75) THEN 
CALL QLii0VR(-5,0) 

GOTO 520 
EM)IF 

IF (L .EQ. 77) THEN 
CALL QLMCM^(5,0) 

GOTO 520 
EM)IF 

IF (L .EQ. 00) THEN 
CALL QLliDVR(0,-5) 

GOTO 520 
EM)IF 

IF (L .EQ. 71) THEN 
CALL GLREAD(XR,YR) 

CALL (lJiDV^(320-XR,17^YR) 

CALL QITDR(XR,YR,XFIN,YFIN) 

X2(l) = XFIN 
Y2(l) = YFIN 
CALL QTPO.(0,N7,X2,Y2) 

CALL QLSTAT(O) 

GOTO 530 
ELSE 
GOTO 510 
EM)IF 
EMHF 

IF (YNATE .EQ. 'VOLTS') THEN 
GOTO 700 
EMIIF 

K6 = 1 
K7 = N3 

DO 600 I = 1,N3 
X = TX(I) 

Y=TY(I) 

IF (((X .GT. XST) .PM}. (X .LT. XFIN)) .AM). ( (Y .GT. 

+ .AM). (Y .LT. YFIN))) ThEN 

K6 = I 
GOTO 605 
ENDIF 
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600 

605 



610 

630 



670 



700 

715 



720 

730 



770 



790 



CONriNUE 
DO 610 I = K6,N3 
X = TX(I) 

Y = TY(I) 

IF ((<X .LT. X5T) .CF. (X .BT. XFIN)) 

+ .CF. (Y .GT. YFIN))) 7FEN 

K7 = I 
GOTO 630 
EM)IF 

continle; 

N5 = K7-K6 
DO 670 I = K6, K7 
YDAT(I-K6<-1) = TY(I) 

XDAT<I-K6+1) = TX(I) 

CUvITINUE 
KO = 1 
K9 = N4 
K2 = N4 
K3 = 1 

DO 700 I = 1,N4 
X = TXl(I) 

Y = 7YKI) 

IF (((X .GT. XST) .PIVO. (X .LT. XFIN)) 
+ .«M). (Y .LT. YFIN))) TVEN 

K9 = I 
GOTO 715 
ENDIF 
CONTINUE 
DO 720 I = K8,N4 
X = TXl(I) 

Y = TYl(I) 

IF (((X .LT. XST) .OR. (X .GT. XFIN)) 

+ .OR. (Y .GT. YFIN))) T>CN 

K9 = I 
GOTO 730 
EM)IF 

cmriNUE 
DO 770 I = KB,K9 
YDATl(I-KB^^l) = TYl(I) 

XDATi(I-K9^1) = TXKI) 

CONTINUE 
N6 = K9 - K8 
GOTO 225 
N5 = O 
N6 = 0 

DO 785 I = 1,N3 
X = TX(I) 

Y=TY( I ) 

IF (((X .GT. XST) .fMi. (X .LT. XFIN)) 



.OR. ((Y .LT. YST) 



.PMy. ((Y .GT. YST) 



OR. ((Y .LT. YST) 



.fiM). ((Y .GT. YST) 
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705 

790 

000 



830 

050 



+ .P#vO. (Y .LT. VFIN))) TVEN 

N5 = N5 +1 
XDAT(N5) = TX(I) 

YDAT(lve) = TY(I) 

ENDIF 
CONTIMJE 
DO 790 I = 1,N4 
X = TXl(I) 

Y=TY1(I) 

IF (((X .GT. XST) .PInO. (X .LT. XFIN)) .PM). ( (Y .GT. 
+ .PM). (Y .LT. YFIN))) 7V04 

N6 = N6 -^1 
XDATKN6) = TXl(I) 

YDATKNt,) = TYl(I) 

EM)IF 
CGNTIMJE 
M = O 
GOTO 225 

OJEB = 'DO YOU WPNT TO PLOT' 

CPLL QGTXT(^,OUE5, 3, 1,335,0) 

QLES = 'THIS GRPFH Y OR N ?' 

CPLL QOTXT(28,QLES,3, 1,320,0) 

CPLL GBEEP 
CPLL QIN<EY(I,M0) 

IF ((MB .EQ. 09) .OR. (MB .EQ. 109)) Th€N 
GOTO 830 

ELSE 

IF ((MB .EQ. 78) .GR. (MB .EQ. IIO) ) THEN 
GOTO 050 
ELSE 

GOTO 0OO 
EM)IF 
ENDIF 

CPLL FLOT(XDAT,YDAT,N5,XDATl,YDATl,N6,XNPIt,YNf#t,TITlJE, 
^-XMAX , XMIN,YmX ,YMIN, XST , XFIN, YST , YFIN, XMAJ0R,YMPJGR,TITLE2, 
4FPCT,NPf-E) 

IF (YNP^ .EQ. 'VOLTS') TFEN 
GOTO 390 
EM)IF 

IF (YNPIt .EQ. 'Dv/Di') TVEN 
GOTO 380 
ENDIF 

DO 8^,0 I = 1,2000 
TX(I)=0 
TY(I)=0 
TXl(I) =0 
TYl(I) = O 
XDAT(I) = O 
XDATl(I) = O 
YDAT(I) = O 
YDATl(I) = O 
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YST) 



860 



CCNTINLE 
COLL QS-CDe(16) 
CflLL IS13DE(3) 
900 FORMAT(Al) 

OSO 
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* 
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* 

t 

t 

t 

t 

t 

t 
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t 

* 
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THIS SLSFOJTINE RJDTS THE GRPRTB AS SHDWM ON THE SCREEN OJRING 
THE CPELING SUBFOJTINE GRf*^l. 

THE mJCRITY CF TTE PRDGRW IS INVOLVED IN SETTING LP PfOEHENTS 
TO BE PASSED TO SLBRDLf^I^E COSTTAirED WITHIN TTE PLOTHATICS 
LIBRWY. THIS MANUAL SHOULD EE CUEULTED IN ORDER TO DETER1INE 
TTE USE OF TTE ARGLEIENTS OTTER TTV^ THOSE DEFIIED BELOW 



XDAT : 

YDAT : 
XDATl : 

YDATl : 
TITLE : 
^ TLE2: 
XNA^E : 
aWLE: 



TTE ARRAY CONTAINING TTE X VALLES OF TTE FIRST CURVE 
TO BE PLDTTED- 
TIE FIRST CURVE'S Y VALUES 

TTE ARRAY CONATINING TTE Y V¥EUES GF TTE SECOND CURVE 
TO BE PLOTTED, 

TTE SECOM) CURVE'S Y VALUES 
TTE TITLE GF TTE PLOT 
AN DISCRIBTION OF TTE GRAPH 
TTE X AXIS TITLE 

TTE Y AXIS TITLE PRCPERLY POSITION AND CATENTATED 



SUBRDUTIIVE PLDT(XDAT,YDAT,^FT1, XDATl, YDATl, NPT2,XNArE,YmE 
+ ,TITLE, XMAX , XMIN, VmX ,YMIN, XST, XFIN, YST, VFIN,XMAJOR,YMAJOR 
+ ,TITLE2,FACT, NAIE) 

REALM XDAT(2000),YDAT(2000),XDAT1(2000),YDAT1(2000) 

INTEGER»2 JCOL( 10) , JROW( 10) ,N^1 ,rsPT2 
INTEGER 13,14, 15, 16, 17, IS, 111, 112, MB 

CHARACTER'S ^¥^,YNA^E♦5,X^IA^E»ll,TITlJEfl5,SAIETE1[12,TITLE2^32 
>, FACTf6,A4*l, QUES<2B 
A4 = ' ' 

SATPLE = YN^:?E//A4//FACT 

XCRG = XST 

YCRG = YST 

DTX = ABS(XMAX-XMIN) 

DX = DTX /MAX(^PT1,^PT2) 

DTY = ^^(YMAX-YMIN) 

IGPT = O 
NKNT = 5 

fvojrs = O 

LABEL = -1 
MEC = 3 
ISIZE = 1 
ISTEP = 2 
IPORT = 1 

CALL ZINIT( ISIZE, ISTEP, IPGRT, BATCH) 

JCL = 20 
JCR = 620 
ICL = 10 
ICR = 320 

13 = 12 

14 = 150 
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15 = 236 

16 = 336 

17 = Zm 

18 = 322 
111 - 386 
112 = 182 

CALL ZCCW(JCL,ia_,16) 

CALL ZCCWV(JCR,ICR,16) 

CALL ZCONV( 13, 14,16) 

CALL ZCO^( 15,16,16) 

CALL ZCONV( 17, 18,16) 

CALL ZCONV( 111, 112,16) 

XMINl =XMIN -.2*DTX 
XMAXl = XmX+ .1»(DTX) 

YMAXl = VWX + -KDTY 
YMINl = YMIN - -l«DTy 
rCEC = 2 

CALL ZPLOT( JCL, JCR, ICL, ICR, XMINl , XMAXl , YMINl , YMAXl , XCFG, YCRG) 
CALL ZFrrOI(X8T,Y5T,Ll,L2) 

CALL ZFIT0I(XF1N,VF1N,L3,L4) 

CALL ZSETIJP(-1, 1,0,1) 

MINOR = 1 

CALL ZXAXIS(XSr,XFIN,XMP^CR, MINOR, LABEL, NDEC) 

CALL ZYAXIS(Y8T,YFIN,YMAJCR,MIN0R,LABEL,NDEC) 

CALL ZXTICS( 10,JCOL, JRDW) 

CfLL ZYTICS(lO,JC0L,JRDW) 

CALL ZLI^E:(L1,L2,L1,L4,1) 

CAiJ_ ZLI^E(L1,L4,L3,L4,1) 

CALL ZLIf€(L3,L4,L3,L2,l) 

CALL ZLI^E(L3,L2,L1,L2,1) 

CALL ZFU 
□UES = ' ' 

CALL QGTXT(28,QUES, 3, 1,335,0) 

CALL GGTXT(28,QUES, 3, 1,320,0) 

CALL QGTXT(28,GLE5, 3, 1,305,0) 

QLES = IF DKIRED INSERT A PEN- 
CALL QGTXT(28,QUES, 3, 1,335,0) 

□UES = 'TTEN PRESS ANY KEY- 
CALL QGTXT(^,aUES, 3, 1,320,0) 

CALL QBEEP 
CALL Q1M<EY(I,MB) 

CALL ZTABL(l,KFri,XDAT,YDAT) 

CALL ZFD 

□UES = - IF DESIRED INSERT A PCW PEN- 
CALL QGTXT(28, CUES, 3, 1,335,0) 

CUES = -TI-EN PRESS ANY KEY- 
CALL QGTXT(28,QLJES,3, 1,320,0) 

CALL rJHFFP 
CALL QIM<EY(I,M8) 

C:ALL ZTAH_(1,M=T2,XDAT1,YDAT1) 

CALL ZFU 
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D5LL QIM<EY(I,MB) 

CUES = • IF DESIRED INSERT A ^CW PEN' 
CiaLI. C3GTXT(28,tUES, 3, 1,335,0) 

CUES = ' TT€N PRESS PNY KEY' 

CALL QGTXT(28,ajES,3, 1,320,0) 

CALL CfflEEP 

CALL ZPTXTA(ll,XNArt,l) 

CALL ZSEJLF(2,1,0,1) 

CALL ZIW( 1,1, 1760, 1400,0) 

CALL ZPTXT(15,TnUE,l,I5,I6) 

CALL ZPTXT(8,NArt,l,Ill,I6) 

CALL ZPTXT(32,TITLE2,1,I5,I8) 

CALL ZEXTXT(0,0,0,0,0,1) 

CALL ZPU 
CALL ZF¥^(I3,I4) 

CALL ZLABEL(12,SAWJE) 

CALL ZFINIS 
END 
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» THIS IS Tl-E FIFST SUBR□UTI^E CPi I FT) W-EN USING TVE CLBIC 

* SPLirvE METHTO. 

« 

* M3ATA 

* XDAT 

« 

t YDAT 

« 

* XDATl 

* YDATl 

* DIV 

« 

* 

$ n-E INTEGER VPRIABLES ACTUPLLY USED IN THIS SLBR□UTI^E AFE 

* USED TO DEFINE TFE REGIGTS WITHIN BCJIH BRF1M>ES AVAILAIiE 

♦ FOR USE 

» 

SUBROUTirE SLXPEl(MDATA,N,J,NfirE,Al,A2,Cl,C2) 

REAL»4 XDAT(2000),YDAT(2000),XDAT1(2000),YDAT1(2000) 

REAL*8 f«ATA(2000,2),DIV(l(XX)) 

INTEGBR*2 I , J ,N,K,f1I ,N3,N4,N1 ,K1 , A1 ,A2,C1 ,C2 
CHfifVCTER NPfE*8,VmE»5,TITIJE2«32,TITLE3»50 
N1 = N 

YNflIE; = 'Dv/Di' 

TITLE2 = TFFEL SLOPES: CUBIC SPLI^E' 

TITLE3=' ' 

C TPE CATHODIC BRfiNCH IS CPECKED FOR AN INCREASE IN LOG 
C CURRENT DENSITY WHICH COULD CAUSE TPE CALCULATION OF A 
C DERIV/ATIVE OF I^FINTE SLOPE 

DO 20 I = J-1,2,-1 

IF (M3ATA( 1-1,2) .UE. rt)ATA(I,2)) GOTO 25 
20 CCNTirsLE 

K = MAX((J-250),1) 

GOTO 26 

25 K = I 

26 MI = J-9 
J1 = K 

N3 = ^®S(K-MI) 

M = 7 

CALL CSPLIN(rt)ATA,N,M,DIV,K,N3) 

DO lOO I = l,N3-3 
XDAT(I) = I«ATA(J-4-I,2) 

YDAT(I) = DIV(N3-I+1) 

100 CONTIMJE 



: TIE ARRAY CCNTAINI^G POTENTIAL AND LOG LLRFENT DE^eITY 
: TIE «EAY CONTAINING TIE LOG CURRENT VALLES FDR TIE 
CATHODIC BRANCH 

; TIE ARRAY CONTAINING TVE DERIVATIVES FOR TIE CATHODIC 
BRANCH 

: S«E A XDAT BUT FOR TIE ANODIC BRANCH 
: SATE AS YDAT BUT FOR TIE ANODIC BRANCH 
: TRE ARRAY CONTAINING TIE DERIVATIVES WHICH IS BEING 
RETURN FROM SUBSEQLENT CALLED SUBROUTI^ES 
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C THE fiMX)IC BRftrCH IS Q-GCKED FDR LUG CURRENT VmJEB WHICM 
C nULD CAUSE TT€ CALJCLLATION GF A CERIVATILE OF IM^INITE SLCPE 

DO 120 I = J+1,N1 

IF (rCATA(I,2) -LE. r^TA( 1-1,2)) GOTO 125 
120 CONTIMJE 

K = MIN((J+250),(N1)) 

GOTO 126 

125 K = I 

126 ai = J 

MI = K - 9 
N4 = ABS(MI-J) 

CALL CSPLIN(M3ATA,N,M,DIV,J,N4) 

DO 200 I = 1,N4 
XDATl(I) = MDATA(J1+I+4,2) 

YDATl(I) = DIV(I) 

200 CGNTINUE 

CALL GRAPHIC XDAT,YDAT,N3,J,NA^E,XDATl,YDATl,N4,VNA^t,TnlJE2, 
+A1 , A2,C1 ,C2,TITL£3) 

EM) 



128 



* THIS SLBR□LrrI^E SETS UP THE MATRIX CF CCFFFICIEWS TO 

* EE SOLVED IN GREER TO DETERMI^E ThE DERIVATIVE LBING THE 

* CUBIC SR_I^E hETHOD, 

« 

* rEIATA : TTE PRRAY CONTAINING POTENTIAL AM) LOB CLFRENT DENSITY 

* CWT : Tl-E 7x7 MATRIX CONTAINING TTE COEFFICIENTS 

* YMAT : Tl-E 7x1 VECTCR CGNTAINING TTE RIGHT HAND SIDE CF 

* TTE SYSTEM CF EGUATICNS. WEN fETURED FRCM TTE 

» SGLVING SUBROLfTINE IT CONTAINS TTE SOLUTION. 

SUBROUTIIE CSPLI^E(M)ATA,J,M,DIV,K,MI) 

REAL*8 XDAT(IOOO), YDAT (1000) ,«)ATA( 2000,2) ,YMAT( 7) 

REAL*8 CMAT(7,7),AI,BI,CI 
REALMS DIV(IOOO) 

IMIEGER»2 N,MI,K,J 
DO lO I = 1,7 
DO 5 J = 1,7 

CMAT(I,J) = O.DO 
5 CONTINUE 

lO (XNT INUE 
M = 7 
K = K-1 

DO 200 N = 1,MI 
20 K = K + 1 

CMAT(1,1) = 2.D0*(M)ATA(K+2,2)-M3ATA(K,2) ) 

01AT(1,2) = M)ATA(K+2,2) - MDATA(K+1,2) 

CMAT(7,7) = 2.DO*(MDATA(K+e,2)-M)ATA(K+6,2) ) 

CMAT(7,6) = M)ATA(K+8,2)-M)ATA(K+7,2) 

YMAT ( 1 ) = ( ( MDATA ( K+2 , 1 ) -M3ATA ( K+ 1 , 1 ) ) / ( M3AT A ( K+2 , 2 ) -MDATA 

+ {K+1,2))) - ((MDATA(K+1,1)-M)ATA(K,1))/(MDATA(K+1,2)-M)ATA 
+(K,2))) 

YMAT(l) = YMAT(1)*6. 

YMAT(7)= ((MDATA(K+8,1)-M)ATA(K+7,1) )/(^f)ATA(K-^e,2)-^l)ATA 
+ (K+7,2))) - ((M)ATA(K+7,l)-M)ATA(K-^A,l))/(M3ATA(K+7,2)-M)ATA 
+(K+6,2))) 

YMAT(7) = YMAT(7)*6. 

DO lOO I = 2,6 
K = K +1 

CMAT(I,I) = 2.D0»((M)ATA(K+2,2)-MDATA(K,2))) 

CMAT( 1,1+1) = M)ATA(K+2,2)-FCATA(K+l,2) 

CMAT( 1,1-1) = M)ATA(K+1,2) - M)ATA(K,2) 

YMAT(I) = ((M)ATA(K+2,1)-M)ATA(K+1,1))/CMAT( 1,1+1)) 

+-( (M)ATA(K+1 , 1 )-M)ATA(K, 1 ) ) /CMAT( I , I-l ) ) 

YMAT(I) = YMAT(I)«6. 

100 CONTirvLE 
K = K - 5 

CALL LIN5Y1(CT1AT,YMAT,M) 

AI = YMAT(5) - YMAT(4) 
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AI = AI/((6-D0)*((riDATA(K^^,2)--I^TA(K+3,2)))) 

B£ = VMAT(4)/2.DO 

Cl = (r^TA(K+4a)~^^TA(K+3,l) )/(^CATA(K+4,2)-^I)ATA(K+3,2) ) 
Cl = CI-(l/6)<(((2tYMAT(4))+YriAT(5))<(MDATA(K+4,2)HiDATA(K+ 
+3,2))) 

DIV(N) = Cl 
200 CONTINUE 
END 



130 



* THIS SL^aRXm^E SOLVES THE 7x7 SYSTEM OF EOUATICNS GET'EFATED 

* BY CEPLIN.FOR. TEH FROGmi WAS CJTIGINALLY FF01 THE TEXTBOOK 

* "FOFnEWN 77 for Engineers and Scientists" BY NVHDFF FND 

* LEESTMA, PAGES 264-265. TFE MAIN MCDIFICATION IS TFE ELIMIEWTION 

* OF TEE VARIAa_E DI^E^BIONI^G OF ARRAYS PRESENT IN TEE ORIGINAL 

$ PROGRAM. 

SLJBROL^■I^E LINBYl (COMAT, YMAT,N) 

INTEGER N, PIVOT 

R0^»8 ALG(7,8), C0MAT(7,7) , YMAT(7) ,X(7) ,MJ_T,TEEP 
DO 20 I = 1, N 
DO 10 J = 1,N 

AOB(I,J) = CCMAT(I,J) 

10 CONTINUE 

20 CONTIMJE 

DO 30 I = 1,N 

AUG(I,N+1) = YMAT(I) 

30 CONTINUE 

DO 70 I = 1,N-1 

IF (AUG(I,I) .EQ. O.) TEEN 
PIVOT = O 
J = I+l 

35 IF ((PIVOT .ED. O) .AND. (J .LE. N) ) TEEN 

IF (AUG(J,I) .^E. O. ) PIVOT = J 
J = J+1 
GOTO 35 
EMDIF 

IF (PIVOT .EQ. O) TEEN 
FRINT», 'SINBULPR* 

STOP 

R 

DO 40 J = 1,N+1 
TEEP = AUG(I,J) 

AUG(I,J) = AUG(PIVDT,J) 

AUG(PIVDT,J) = TEEP 
40 CONTINUE 

ENDIF 
ENDIF 

DO 60 J = I+1,N 

MJLT = -AUG(J,I)/AUG(I,I) 

DO 50 K = I,N+1 

AUG(d,K) = ALJB(J,K) + MJLT««JG( I ,K) 



50 


CONTIMJE 


60 


CONTINUE 


70 


CCiNTINUE 




X(N) = ALG(N,N+1)/AUG(N,N) 




DO 90 J = N-1,1,-1 



X(J) = AUG(J,N4-1) 
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00 

90 

100 



IX) 00 K = Jf-1,N 

X(J) = X(J)-(PUG(J,K)-*X(K) ) 
OlSfriNUE 

X(J) = X(J)/AUG(J,J) 

(XTsfriMJE 
CO lOO I = 1,N 
YMAT(I) = X(I) 

CONTINUE 

EM) 
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* THIS aJERCUriNE COv^EFfTS A MJTERIC VfiRIBLE (X6) CF FORMAT 

* +/- XX.XXXX TO A OWVCTER VWIBUE WHOSE WIDTH IS 8 CHPRACTERS. 

» THIS IS ^ECESSARY AS TFE R_DmATICS LIBRARY IS NOT EOUIFFED WITH 

* A SLJBR□Lm^E FOR PLOTTING MJTERIC VflRIBLK. 



SLBRCUTINE CCNV(A7,X6) 

INTEGER 01,02, 03,04,05, 06 

CHARACTER*! A6(0:9), Cl ,C2,C3,C4,C5,C6,C7,C8,A7*8 
A6(0) = *0' 

A6<1) = 1' 

A6(2) = '2' 

A6(3) = '3' 

A6(4) = '4' 

A6(5) = '5' 

A6(6) = 'h' 

A6(7) = ’7‘ 

A6(8) = '8' 

A6(9) = '9' 

Cl = ' ' 

C2 = ' ' 

C3 = ' 

C4 = ' 

C5 = ' ' 

C6 = ' 

C7 = ' ‘ 

IF (X6 .LT. 0.0) Cl = 

X6 = AB5(X6) 

XI = X6/10. 

01 = IFIX(Xl) 

X2 = X6-(10.»01) 

02 = IFIX(X2) 

X3 = X6 - (10.*01) - 02 
X3 = X3/.1 

03 = IFIX(X3) 

X4 = (X3-03)/.l + .OOl 

04 = IFIX(X4) 

X5 = (X4-04)/.l 

05 = IFIX(X5) 

X7 = (X5-05)/.l 

06 = IFIX(X7) 



10 


I = 


0,9 










IF 


(01 


.ED. 


I) 


02 = 


A6( 


I) 


IF 


(02 


.ED. 


I) 


C3 = 


A6( 


I) 


IF 


(03 


.ED. 


I) 


C4 = 


A6( 


I) 


IF 


(04 


.ED. 


I) 


C5 = 


A6( 


I) 


IF 


(05 


.ED. 


I) 


Ch = 


A6( 


I) 
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10 



IF (06 .EQ. I) C7 = A6(I) 
casrriNLE 
C8 = ' . ' 

A7 = C1//C2//C3//C8//C4//C5//C6//C7 

RETURN 

END 
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* THIS SLSTOLfTIfvC; SETS SCTE CF Tl-E CHWACTER STRINGS LEED 

* lO-EN USING Tl-E CENTRAL DIFFEFENCE 1‘ETKDD 

SUERDUTIfvE SLXFE(MDATA,N,J,NP<«E,A1,A2,C1,C2,N1,N2,N3,N4) 
FEPL*4 XDAT ( 2000 ),YDAT( 2000 ),XDAT1( 2000 ),YDAT1( 2000) 

REAL*8 ^fl3ATA(2000,2) 

INTEGER»2 I ,J,N,A1,A2,C1,C2 

CHARACTER ^W^E»8,YNA^E»5,TITLE2»32, LIST*50,TITLE3*50 
YNPME = Dv/Di' 

TIT1_E2 = TAFEL SLOPES: CENTRAL DIFFERENCE' 

TITLE3=' ' 

CALL DATAiEL ( f«ATA , J , YDAT , YDATl , XDAT , XDATl , N2 , N1 ) 

CALL 6FVRHK XDAT, YDAT,N3,J,!v*^, XDATl, YDATl ,N4,YT^=<*E,TITL£2, 
+A1,A2,C1,C2,TITLE3) 

END 
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t THIS SUBROUTINE CONDUCTS NUnEJRICPi^ DIFFERENTIATICN USING 
♦ THE 4 POINT CENim. DIFFEPNCE rETHDD, 



t 

t 

« 

t 

t 

♦ 

* 

t 

t 

« 

m 

* 

t 

t 

t 



I^tA = PJRRAY CONTAINING POTENTIPiL AfvO LOB I DETJSITY 
SLOPE = ^«RAY CONTAINING DE/D(LDGI ) AND LOG I 
FSLCFE = CENTRAL DIFFEPNCE DIFFERENTIATION 
ECORRl = UPPEHLIMIT OF ANODIC/CATVODIC DISTRIBUTION 
EC0FR2 = LOUEFI.IMIT OF ^VOIIC/CATHDDIC DISTRIBUTION 
N = NUMBER OF DATA POINTS 

J = LOCATION WITHIN MDATA OF ECORR 

I = CATHODIC BRANCH LOOP COUNTER 

K = ANODIC BRANCH LOOP COUNTER 

H = DIFFE3RENTIAL VALUE IN POTENTIAL 

II = VALUE OF LOG I AT E-2H 

12 = VALUE OF LOG I AT EH-< 

15 = VALUE OF LOG I AT E 

14 = VALUE OF LOG I AT E+H 

15 = VALUE OF LOG I AT E+2H 

DEDI = SLOPE VALUE OF 13 



SUBROUTINE DATADEL (MDATA, J ,YDAT, YDATI,XDAT,XDAT1,N2,N1 ) 
REALfS MDATA(200O,2),Xl,X2,X3,X4,X5,H,FSL0PE,DEDI 
FEAL»4, XDAT1(1000),YDAT1(1000), YDAT( 1000) ,XDAT( lOOO) 
INTEGERt2 J,K,I 

FSLCFE (H,X1,X2,X4,X5)=(XH-(8.D0<X4)-(8.D0»X2)-X5)/(12,D0»H) 
C PERFDR1ING CALCULCATIONS ON CATHODIC BRANCH 



K = O 

DO lOO I = J-2,Nl+2,-l 
K = K+1 

H = MIATA(I,1) - MDATA(I-1,1) 

XI = M}ATA( 1-2,2) 

X2 = MDATA( 1-1,2) 

X3 = MDATA(I,2) 

X4 = MDATA( 1+1,2) 

X5 = MDATA( 1+2,2) 

DEDI = l.D0/FSLQFE(H,Xl,X2,X4,X5) 
YDAT(K) = DEDI 
XDAT(K) = MDATA(I,2) 

100 CONTINUE 
I = O 



C PERF0R1ING CALCULATIONS ON TVE ANODIC BRANCH 



200 DO 300 K = J+2, N2-2 
1 = 1 + 1 

H = riDATA(K+l,l) - MDATA(K,1) 
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XI = MDATA(K-2,2) 

X2 = MDATA(K-1,2) 

X3 = t*1DATA(K,2) 

X4 = MDATA(K+1,2) 

X5 = ^CATA(K+2,2) 

DEDI = l.DO/FSLOPE(H,Xl,X2,X4,X5) 
YDATl(I) = DEDI 
XDATl(I) = MDATA(K,2) 

300 CONTINLE 
400 RETUFaM 
EM> 
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